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FOREWORD 


This  report  represents  the  results  of  the  efforts  expended  in 
performance  of  Contract  F33615-69-C-1L2L t  ’’Development  of 
Stability  and  Control  Prediction  Methods  for  Stoppable  Rotor 
Aircraft.”  The  work  was  performed  by  Bell  Helicopter  Company 
under  Project  No.  8219.  It  was  sponsored  by  the  Air  Force 
Flight  Dynamics  Laboratory,  Air  Force  Systems  Command,  from 
December  1968  through  February  1970.  Mr.  Charles  L.  Livingston 
was  the  Bell  Helicopter  Company  Project  Engineer.  Mr.  Robert 
Nicholson  was  the  Air  Force  Project  Engineer. 

This  final  report  is  presented  in  four  volumes.  The  first 
describes  the  mathematical  model  and  the  methods  used  to 
calculate  stability  characteristics.  They  are  of  sufficient 
complexity  that  a  digital  computer  is  necessary  for  the  solu¬ 
tion  of  the  equations.  Hie  second  volume  presents  the  re¬ 
sults  of  sample  computations  and  discusses  input  and  output 
formats  and  good  user  techniques.  The  third  volume  describes 
the  computer  program  while  the  fourth  volume  contains  Appen¬ 
dices  which  are  computer  generated  documentation  of  the  pro¬ 
gram. 

Hie  authors  gratefully  acknowledge  the  assistance  of  Messrs. 

B.  L.  Blankenship  and  G.  L.  Livingston  of  the  Bell  Helicopter 
Company  Flight  Technology  Section  in  the  development  of  the 
mathematical  model. 

"his  technical  report  has  been  reviewed  and  is  approved. 


Chief,  Control  Criteria  Branch 

Flight  Control  Division 

Air  Force  Flight  Dynamics  Laboratory 


ABSTRACT 


This  volume  presents  all  documentation  available  to  aid  the 
user  of  the  computer  program  developed  in  this  work.  The 
input  format  section  provides  an  explanation  of  all  of  the 
quantities  input  to  the  computer  program.  Many  of  the  inputs 
are  defined  by  equations  showing  how  they  function  in  the  pro¬ 
gram.  This  makes  the  use  of  the  inputs  as  clear  as  possible. 
Four  typical  sets  of  input  data  are  included  as  working  ex¬ 
amples  . 

The  output  guide  gives  a  thorough  discussion  of  all  of  the 
forms  of  computer  output  obtained  by  the  user. 
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SECTION  I 


INTRODUCTION 


The  purpose  of  this  document  is  to  inform  the  reader  of  the 
various  capabilities  of  computer  program  ASAJOL  and  to  pro¬ 
vide  the  information  necessary  for  the  assembling  of  data 
into  the  proper  form  and  format  for  a  successful  execution  of 
the  program. 

The  first  section,  the  Input  Format,  lists  the  input  data  in 
a  sequence  indicating  the  position  of  each  parameter  in  the 
data  deck.  Where  possible,  a  one  line  description,  including 
units,  is  used.  Some  of  the  inputs  cannot  be  so  described. 

In  these  cases,  a  FORTRAN  symbol  or  some  other  symbol  is  used. 
Where  a  FORTRAN  symbol  is  used,  the  following  conventions  are 
observed : 

1.  If  the  first  letter  is  I,  J,  K,  L,  M,  or  N,  the  number 
input  should  end  in  the  rightmost  column  of  the  field 
specified  and  should  not  contain  a  decimal  point;  i.e., 

"I"  format  code. 

2.  If  the  first  letter  is  not  I,  J,  K,  L,  M,  or  N,  the 
number  input  should  begin  in  the  leftmost  column  of  the 
field  specified  and  should  contain  a  decimal  point;  i.e., 
"F"  format  code. 

There  are  five  aerodynamic  surfaces  to  be  defined  by  inputs: 
two  rotors,  wing,  elevator  and  fin.  Since  these  inputs  have 
a  uniform  format  and  usage,  only  one  item  by  item  enumeration 
and  explanation,  under  the  title  Aerodynamic  Inputs,  is 
necessary.  The  appropriate  positions  for  these  five  sets  of 
inputs  are  then  indicated. 

The  first  section,  the  Input  Format,  is  designed  to  be  the 
only  documentation  needed  by  the  very  experienced  user.  The 
less  experienced  the  user,  the  more  inadequate  are  some  of  the 
one  line  descriptions  and  symbols.  The  second  section, 
arranged  in  the  same  order  as  the  first,  should  be  referred 
to  in  such  cases  for  a  more  complete  explanation. 

Ir»  this  document  and  in  the  output,  there  are  references  to 
the  rotors  by  various  names.  Where  practical,  the  names  used 
are  appropriate  to  the  configuration.  The  following  two 
groups  of  names  are  the  ones  used  and  the  names  within  a  group 
may  be  considered  synonymous  with  context  determining  the 
appropriate  word: 


1.  Main,  forward,  right,  first 

2.  Tail,  aft,  left,  second 

Each  general  operation  which  the  program  can  perform  is  illus¬ 
trated  in  the  flow  charts  in  Figures  1  through  10,  both  alone, 
if  possible,  and  in  all  permissible  combinations  with  other 
operations.  Each  operation  or  combination  of  operations  is 
controlled  by  data  supplied  by  the  user.  Thus  each  flow  chart 
begins  with  the  operation  "read  data  deck".  Since  the  amount 
of  data  to  be  read  depends  on  the  operation  or  operations  de¬ 
sired,  a  data  deck  in  this  context  consists  of  a  card  ("NPART 
Card”  or  CARO  01  in  the  Input  Format)  telling  the  program 
which  operation  to  perform  followed  by  the  additional  data 
necessary  to  the  performance  of  that  operation.  In  some 
cases,  the  additional  data  are  contained  on  80  or  more  addi¬ 
tional  cards  (trim,  stability  analysis,  maneuver);  in  other 
cases,  no  additional  data  are  required.  Thus,  a  data  deck 
consists  of  an  "NPART  Card"  and  all  subsequent  cards  up  to, 
but  not  including,  the  next  "NPART  Card".  Obviously,  a  deck 
of  cards  to  be  used  in  one  run  on  the  computer  may  consist  of 
any  number  of  data  decks.  There  are  some  restrictions  on  the 
order  in  which  data  decks  may  be  used  in  a  run;  these  will  be 
pointed  out  in  the  appropriate  places. 

The  second  step  on  a  number  of  the  flow  charts  is  "calculate 
problem  constants".  In  each  operation  containing  this  step, 
a  number  of  quantities,  which  must  be  defined  by  the  input 
data,  remain  constant  throughout  the  performance  of  the  opera- 
tion(s).  This  occurs  because  of  the  desirability  of  having  a 
minimum  amount  of  input  data  and  at  the  same  time  having  maxi¬ 
mum  flexibility  in  program  capabilities. 

In  finding  the  trim  point,  the  program  manipulates  the  pilot’s 
controls  and  fuselage  and  rotor  attitude  in  spac^  *-o  reach 
desired  values  of  the  rotor  flapping  moments  ar  tl  e  resultant 
forces  and  moments  on  the  center  of  gravity. 

Sometimes  these  desired  values  are  all  zero.  In  this  case, 
the  trim  point  is  such  that  the  rotorcraft  would  continue  to 
fly  indefinitely  in  a  straight  line  in  space.  That  is,  main¬ 
taining  an  input  rate  of  climb  or  descent  as  well  as  forward 
and  lateral  velocity. 

Sometimes  these  desired  values  are  not  all  zero.  In  this 
case,  the  use  of  the  word  "desired"  rather  than  "predetermined" 
is  significant  because  the  desired  values  depend  on  the  trim 
point.  Two  different  results  are  possible  in  this  case.  In 
one,  the  rotorcraft  is  in  a  push-over  or  pull-up  condition  at 
a  predetermined  g-level.  This  condition  cannot  be  maintained 


indefinitely  and  should  not  be  attempted  in  conjunction  with 
a  maneuver.  In  the  othe the  rotorcraft  is  in  a  banked  turn, 
either  level  or  spiral,  at  a  predetermined  g-level,  bank 
angle,  or  turn  radius.  This  condition  may  be  used  as  the 
starting  point  for  a  maneuver. 

The  stability  operation  uses  a  trim  point  as  its  initial  con¬ 
dition.  Calculated  from  this  condition  are  the  stability 
derivatives  of  the  rotorcraft,  both  longitudinally  and  later¬ 
ally  in  a  free  flight  condition  with  no  inputs  and  then  with 
a  step  input  from  each  of  the  controls.  Transfer  functions 
are  also  computed. 

A  parameter  sweep  may  be  used  when  either  finding  a  trim  point 
or  doing  a  stability  analysis.  The  most  frequent  uses  are 
for  sweeping  speed  or  ^-level,  but  any  input  or  combination 
of  inputs  used  in  finding  a  trim  point  or  doing  a  stability 
analysis  may  be  swept.  The  assumption  is  made  that  each  data 
point  bears  some  relationship  to  the  previous  one,  and  further 
that  the  trim  point  of  each  one  is  a  good  starting  point  for 
finding  the  trim  point  of  the  next  one.  For  example,  in  a 
speed  sweep,  a  change  of  20  or  30  knots  is  the  most  that  should 
normally  be  used  above  40  knots.  Below  40  knots,  the  maximum 
change  should  not  exceed  10  knots. 

A  maneuver  also  requires  a  trim  point  as  a  preliminary.  Here 
the  trim  point  is  used  to  supply  the  initial  conditions  to  a 
system  of  differential  equations  which  describe  the  behavior 
of  a  rotorcraft  in  a  maneuver  under  the  specified  conditions. 
Various  external  inputs,  or  forcing  functions,  may  be  applied. 
Examples  are  control  movements,  gusts,  and  throttle  chops. 

The  fuselage  and  rotor  blades  are,  in  general,  assumed  to  be 
rigid . 

A  maneuver  with  stability  analysis  is  a  combination  of  two 
operations,  maneuver  and  stability  analysis.  At  times  deter¬ 
mined  by  input  data,  the  maneuver  is  suspended  while  a  sta¬ 
bility  analysis  is  performed.  The  maneuver  is  then  resumed  as 
if  no  interruption  had  occurred  and  continues  until  it  reaches 
the  next  time  point  to  do  a  stability  analysis. 

The  value  of  a  large  number  of  variables  is  preserved  at  each 
time  point  in  a  maneuver.  At  the  conclusion  of  the  maneuver, 
the  user  may  select  any  of  these  variables  to  be  plotted  on  a 
scale  of  his  choice. 

At  the  conclusion  of  a  maneuver  all  of  the  variables  which 
have  been  preserved  at  each  time  point  during  the  maneuver  may 
then  be  transferred  to  a  tape  which  can  be  stored  in  the  tape 
library.  Later  the  data  may  be  recovered  for  any  use  desired. 
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At  the  conclusion  of  a  maneuver,  usually  one  in  which  one  of 
the  controls  or  the  mast  tilt  angle  has  been  varied  sinusoid¬ 
ally,  any  of  the  variables  which  have  been  stored  during  the 
maneuver  may  be  selected  for  analysis.  The  analysis  consists 
of  fitting,  by  least  squares,  a  curve  of  the  form 


F.(t)  =  A.sin(ci)t  +  0.)  +  B. 

11  11 

to  the  data  stored.  Then  any  amplitude  ratios,  Ai/Aj,  and 
phase  angle  differences,  0i  -  0j,  may  be  computed.  Then 
linear  combinations  of  the  variables,  in  the  form 


Fi(t)  =  C±  F..(t)  +  Di  Fk(t)  +  Eit 


may  be  derived. 

Some  sample  data  setups  are  included  to  help  the  novice  user 
get  started  and  to  give  some  insight  into  orders  of  magnitude 
for  some  of  the  inputs. 

The  output  from  the  computer  program  can  be  separated  into 
useful  results  for  the  user  and  messages  resulting  from  an 
error  condition.  Each  type  of  output  has  been  explained  in 
detail.  Some  error  conditions  cannot  be  diagnosed  in  the 
program.  An  example  would  be  an  unrealistic  flight  condition 
which  cannot  be  trimmed.  This  type  of  problem  depends  on  the 
knowledge  and  insight  of  the  user  for  detection  and  correction. 
Some  hints  have  been  given  to  help  the  new  user  in  these 
areas. 
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SECTION  II 
INPUT  DATA 


A.  INPUT  FORMAT 

1  *  Identification  and  Flow  Control  Group 


CARD  00 


CARD  01 


Message  instructing  computing  control  section  about 
alphameric11  °f  output  and  data  deck.  Columns  1-68, 


Col.  1-2 

NPART  (Permissible  v 

1,2,3,7,8,10. 

Col.  4-6 

NPRINT 

Col.  10 

NSCALE 

Col.  13  -  15 

NVARA 

Col.  21  -  25 

AL(l) 

Col.  26  -  30 

AH(1) 

Col.  33  -  35 

NVARB 

Col.  41  -  45 

AL(2) 

Col.  46  -  50 

AH(2) 

Col.  53  -  55 

NVARC 

Col.  61  -  65 

AL(  3) 

Col.  66  -  70 

AH(3) 

CARD 

02 

Col.  3-10 

IPSN 

Col.  11  -  70 

Identifying  Comments 

CARD 

03 

Col.  1  -  68 

Identifying  Comments 

CARD 

04 

Col.  1  -  68 

Identifying  Comments 

CARDS 

11-201  have  the  format  7F10.0 

( 

h 
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2.  Aerodynamic  Inputs 


CARD  A 

YXX*  (1) 
(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

CARD  B 

(8) 
(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

CARD  C 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 


Drag  divergence  Mach  number  for  a=0 

Mach  number  for  lower  boundary  of  supersonic 

region 

Maximum  Cl,  normal  flow,  M  (Mach  number) =  0 

Coefficients  of  Mach  number  in  maximum 
Cl  equation,  normal  flow 

Maximum  Cl,  reversed  flow,  M=0 


I  Coefficients  of  Mach  number  in  maximum 
Cl  equation,  reversed  flow 

Tail  boom  bending  coefficient  (/lb) 

Cd  for  a =0,  M=0 

|  Coefficients  of  a  in  non-divergent  (/deg) 
J  drag  equation  (/deg2) 


Coefficient  in  supersonic  drag 
eauation 

Maximum  non- divergent  Cq 

IVo-dimensional  lift  curve  slope  (/deg) 

(if  zero,  data  tables  are  used) 

Aspect  ratio 


*YXX  stands  for  YMR,  YTR,  YWG,  YEL,  or  YFN  with  context 
determining  which  one.  Card  numbers  A,  B,  and  C  are  used 
here  only  to  indicate  order.  Actual  numbers  vary  from  group 
to  group  and  are  given  in  the  proper  place  in  each  group. 
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3.  Fuselage  Group 


CARD  11 
XFS 

(1) 

Gross  weight 

(2) 

Station  line 

(3) 

Buttline 

(4) 

Waterline 

(5) 

Station  line 

(6) 

Buttline 

(7) 

Waterline 

Location  of  fuselage 
aerodynamic  center 

Location  of  center 
of  gravity 


(lb) 

(in) 

(in) 

(in) 

(in) 

(in) 

(in) 


CARD  12 

(8)  Fuselage  rolling  inertia,  Ix 

(9)  Fuselage  pitching  inertia,  Iy 

(10)  Fuselage  yawing  inertia,  I2 

(11)  Fuselage  product  of  inertia,  Ixz 

(12) 

(13) 

(14) 


(slug-ft^) 
(slug- ft2) 

(8lUg-ft2) 

(8lug-ft2) 


CARD  13 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 


Coefficients  in 
moment  equation 
Coefficients  in 
moment  equation 


fuselage  pitching  (ft^) 

(ft3/deg) 

fuselage  yawing  (ft3) 

(ft3/deg) 


Coefficients  in  fuselage  lift 
equation 


(ft2) 

(ft2/deg) 


CARD  14 

(22) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


Coefficients  in  fuselage  drag 
equation 

Coefficients  in  fuselage  side 
force  equation 


(ft2) 


( ft2/ deg) 
(f t2/deg2) 

(ft2/deg2) 

(ft2) 
(f t2/deg) 
(f t2/deg2) 


CARD  15 

(29)  Pylon  weight  (lb) 

(30)  Station  line  |  Location  of  pylon  center  (in) 

(31)  Waterline  (  of  gravity  (in) 

(32)  Pylon  differential  flat  plate  drag  area  (ft2) 

(33)  Distance  from  mast  pivot  point  to  pylon 

aerodynamic  center  (ft) 

(34) 

(35) 
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4.  Main  Rotor  Grout 


CARD  21 
XMR 


Number  of  blades 

Number  of  radial  aerodynamic  stations 


Radius 

Chord 

Coning  spring  constant 
Flapping  stop 


(ft) 

(in) 

(ft-lb/deg) 

(deg) 


CARD  22 


(8)  Station  T  .  .  c  .  . 

(9)  Buttline  Location  of  shaft  pivot 

(10)  Waterline  point 

(11)  Mast  tilt  angle 

(12)  Mast  length 

(13)  Rotor  to  engine  gear  ratio 

(14)  Fold  rotor  C.G.  shift  switch 


(in) 

(in) 

(in) 

(deg) 

(ft) 


CARD  23 


CARD  24 


(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 


(22) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


Virtual  hinge  point 
Flapping  hinge  offset 
F/A  flapping  spring  rate 
Lateral  flapping  spring  rate 
Spring  rate  for  focused  pylon 

Hub  extent 


Precone 

Pitch- cone  coupling  ratio 
Pitch-flap  coupling  ratio,  83 
Drag  coefficient  for  hub 
Radial  flow  factor 
Coefficient  for  tip  vortex  effect 


(ft) 

(ft) 

(ft-lb/deg) 

(ft-lb/deg) 

(deg/lb) 

(ft) 


(deg) 


(deg) 


CARD  25 

26  Twist 

27 


CARD  28 
29 
2A 

CARD  2B 
2C 
2D 


Blade  weight  distribution 


YMR(l)  -  YMR( 21) 


(deg) 


(1 bs/iny 


5,  Tail  Rotor  Group 


CARD  31 
XTR 


(1)  Number  of  blades 

(2)  Number  of  radial  aerodynamic  stations 

(3) 

(4)  Radius  (ft) 

(5)  Chord  (in) 

(6)  Coning  spring  constant  (ft-lb/deg) 

(7)  Flapping  stop  (deg) 


CARD  32 


(9)  Buttline  |  Lo?ation  of  shaft  Pivot 

(10)  Waterline  )  Polnt 

(11)  Mast  tilt  angle 

(12)  Mast  length 

(13)  Rotor  to  engine  gear  ratio 

(14)  Tail  rotor  indicator 


(in) 

(in) 

(in) 

(deg) 

(ft) 


CARD  33 

(15)  Virtual  hinge  point 

(16)  Flapping  hinge  offset 

(17)  F/A  flapping  spring  rate 

(18)  Lateral  flapping  spring  rate 

(19)  Spring  rate  for  focused  pylon 

(20) 

(21)  Hub  extent 


(ft) 

(ft) 

(ft-lb/deg) 

(ft-lb/deg) 

(deg/lb) 

(ft) 


CARD  34 

(22)  Precone  (deg) 

(23)  Pitch-cone  coupling  rat ib 

(24)  Pitch- flap  coupling  ratio,  §3 

(25)  Drag  coefficient  for  hub 

(26)  Radial  tlow  factor  (deg) 

(27)  Coefficient  for  tip  vortex  effect 

(28)  Sidewash  coefficient 


CARD  35  ) 

36  [  Twist 

37  ) 


CARD  38  ) 

39  J  Blade  weight  distribution  (lbs/in) 

3A  ) 


CARD  3B  ) 

3C  YTR(l)  -  YTR( 21) 
3D  ) 
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6 .  Wing  Groui 


CARD  41 

XWG  (1) 


CARD  42 


CARO  43 


(6) 

(7) 


Wing  area 

Station  line  T  .  .  ~  .  - 

Buttline  Location  of  center  of 

Waterline  1  Pressur®  of  rl8ht  w«>g 
Geometric  angle  of  incidence  relative 
to  fuselage  centerline 


(ft2) 

(in) 

(in) 

(in) 

(deg) 


(10) 

(11) 

(12) 

(13) 

(14) 


(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 


Main  rotor  induced  velocity  factor 
Coefficient  in  wing  wake  deflection 
equation 

Coefficient  in  >?q  equation 


(deg) 


Coefficient  of  sideslip  in  roll  moment  equation 
Coefficient  of  sideslip  and  Cl  in  roll  moment 
equation 

Coefficient  of  yaw  rate  and  Cl  in  roll  moment 
equation 


Coefficient  of  roll  rate  in  roll  moment 
equation 

Coefficient  of  sideslip  in  yaw  moment  equation 
Coefficient  of  sideslip  and  Cl2  in  yaw  moment 
equation 

Coefficient  of  yaw  rate  and  Cl2  in  yaw  moment 
equation 

Coefficient  of  yaw  rate  and  dC^/d a  in  yaw 
moment  equetion 

Coefficient  of  roll  rate  and  Cl  in  yaw  moment 
equation 

Coefficient  of  roll  rate  and  Cp  in  yaw  moment 
equation 


CARD  44  1 

45  f  YWG(l)  -  YWG(21) 

46  ) 


' 


- 


7 .  Elevator  Grou 


CARD  51 

XEL  (1)  Area 

(2)  Station  line 

(3)  Buttline 

(4)  Waterline 


Location  of  center  of 
pressure 


(5)  Geometric  angle  of  incidence  relative 
to  fuselage  centerline 

(6) 


(7) 


(ft2) 

(in) 

(in) 

(in) 

(deg) 


CARD  52 

(8)  Induced  velocity  factor 

(9)  Velocity  at  which  elevator  starts  to 

enter  main  rotor  wake  (kts) 

(10)  Velocity  at  which  elevator  is  completely 

into  main  rotor  wake  (kts) 

(11)  Coefficient  for  change  of  angle  of 

attack  due  to  wing  wake  (deg) 

(12) 

(13) 

(14) 


CARD  53  ) 

54  YEL(l)  -  YEL( 21) 

55  ) 
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3.  Fin/Rudder  Grou 


CARD  61 

XFN  (1)  Area 

(2)  Station  Line 

(3)  Buttline 

(4)  Waterline 


Location  of  center  of 
pressure 


(5)  Geometric  angle  of  incidence  relative 
to  fuselage  centerline 

(6)  Induced  velocity  factor 

(7)  Sidewash  coefficient 


CARD  62  } 

63  YFN(l)  -  YFNC21) 

64 


9 .  Jet  Groui 


CARD  71 

XJET  (1)  Number  of  jets 

(2)  Thrust  of  right  jet 

(3)  Thrust  of  left  jet 

( 4)  Station  line 

(5)  Buttline 

(6)  Waterline 

(7) 


Location  of  right  jet 
thrust 


(lb) 

(lb) 

(in) 

(in) 

(in) 

(in) 


(8)  Yaw  angle,  body  to  right  jet 

(9)  Pitch  angle,  body  to  right  jet 
1  °  U 


CARD  72 


(deg) 

(deg) 


I 


10.  Bobweight  Group 


CARD  81 

XBW 

(1) 

V 

(2) 

Kx 

(3) 

cx 

(4) 

II 

(5) 

12 

(6) 

13 

(7) 

g  preload 

CARD  82 

(8) 

Mi 

(9) 

mb 

(10) 

m3 

(11) 

M4 

(12) 

(13) 

(14) 

CARD  83 

(15) 

rl 

(16) 

r3 

(17) 

r4 

(18) 

r5 

(19) 

^6 

(20) 

r7 

(21) 

r8 

o 

(deg/sec  ) 
(lb/in) 
( lb-sec/in) 
(lb-in^) 
(lb-in2) 
(lb-in2) 
(g's) 


(lb) 

(lb) 

(lb) 

(lb) 


(  in) 
(in) 
(in) 
(in) 
(in) 
(in) 
(in) 
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* 

,A  * 


Controls  Groui 


a.  Collective 


CARD  91 
XCON 


CARD  92 


(10) 


(11) 


(12) 


(13) 

(14) 


F/A  Cyclic 


Range  of  collective  stick  (in) 

!  Coefficients  for  lower  limit  of 

collective  pitch  as  a  function  (/deg) 
of  An  (/deg2) 

)  Coefficients  for  range  of  col-  (deg) 

)  lective  pitch  as  a  function  of  /9m  (deg) 

Collective  lock  indicator  for  rotor  1 


(/deg) 

(/deg2) 

(deg) 


Root  collective  pitch  of  rotor  1  if 
XCON(7)/0  (deg) 

Coefficient  for  linking  aft/left  rotor 
collective  pitch  to  forward/right  rotor 
collective  pitch  (deg/deg) 

Coefficient  for  linking  stick  posi¬ 
tion  to  geometric  angle  of  incidence 
of  wing  (deg/in) 

Coefficient  for  linking  stick  posi¬ 
tion  to  geometric  angle  of  incidence 
of  elevator  (deg/in) 

Coefficient  for  linking  change  of 
stick  position  to  change  of  auxiliary 
thrust  (lb/in) 

Maximum  rate  for  prop-rotor  collective 
governor  (deg/sec) 

RPM  dead  band  for  prop- rotor  col-  (rpm/rpm) 
lective  governor 


(deg/in) 


(deg/in) 


(lb/in) 


CARD  93 


CARD  94 


(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 


(22) 


Range  of  F/A  cyclic  stick 

F/A  cyclic  pitch  on  swashplate  with 

stick  full  aft 

F/A  range  of  swashplate 

F/A  cyclic  pitch  lock  indicator 

F/A  cyclic  pitch  of  rotor  1  if 

XCON(18)/0 


(in) 

(deg) 

(deg) 


(deg) 


Coefficient  for  linking  forward  rotor 
collective  pitch  to  forward  rotor  F/A 
cyclic  pitch  ( 


(deg/deg) 


L5 


(23)  Coefficient  for  linking  aft  rotor 
collective  pitch  to  forward  rotor 

F/A  cyclic  pitch  (deg/deg) 

(24)  Coefficient  for  linking  aft/left 

rotor  F/A  cyclic  pitch  to  forward/right 
rotor  F/A  cyclic  pitch  (deg/deg) 

(25)  Coefficient  for  linking  right  rotor 
lateral  cyclic  pitch  to  right  rotor 

F/A  cyclic  pitch  (deg/deg) 

(26)  Coefficient  for  linking  left  rotor 
lateral  cyclic  pitch  to  right  rotor 

F/A  cyclic  pitch  (deg/deg) 

(27)  )  Coefficients  for  linking  geo-  (deg/in) 

(28)  |  metric  angle  of  incidence  of  (deg/in2) 
)  elevator  to  stick  position. 


Lateral  Cyclic 


CARD  95 

(29)  Range  of  lateral  cyclic  stick 

(30)  Lateral  cyclic  pitch  on  swashplate 
with  stick  full  left 

(31)  Lateral  range  of  swashplate 

(32)  Lateral  cyclic  pitch  loc  indicator 

(33)  Lateral  cyclic  pitch  of  rotor  1  if 
XC0N( 32)/0 

(34) 

(35) 


(in) 

(deg) 

(deg) 

(deg) 


CARD  96 

(36)  Coefficient  for  linking  right  rotor 
collective  pitch  to  right  rotor 

lateral  cyclic  pitch  (deg/deg) 

(37)  Coefficient  for  linking  left  rotor 
collective  pitch  to  right  rotor 

lateral  cyclic  pitch  (deg/deg) 

(38)  Coefficient  for  linking  right  rotor 
F/A  cyclic  pitch  to  right  rotor 

lateral  cyclic  pitch  (deg/deg) 

(39)  Coefficient  for  linking  left  rotor 
F/A  cyclic  pitch  to  right  rotor 

lateral  cyclic  pitch  (deg/deg) 

(40)  Coefficient  for  linking  aft/left 
rotor  lateral  cyclic  pitch  to 

fwd/right  rotor  lateral  cyclic  pitch  (deg/deg) 

(41)  Lateral  cyclic  stick  -  aileron  coupling 

factor  (deg/in) 

(42) 
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d.  Pedal  and  Mast  Tilt 


CARD  97 

(43)  Range  of  pedal 

(44)  Tail  rotor  collective  or  rudder  angle 
with  pedal  full  right 

(45)  Range  of  tail  rotor  collective  or 
rudder 

(46)  Rotor  2  collective  lock  indicator 

(47)  Rotor  2  root  collective  pitch  if 
XCON(46)/0 

(48)  Rudder  lock  indicator 

(49) 

CARD  98 

(50)  Coefficient  for  linking  right  rotor 
collective  pitch  to  rudder  position  (deg/deg) 

(51)  Coefficient  for  linking  left  rotor 
collective  pitch  to  rudder  position  (deg/deg) 

(52)  Coefficient  for  linking  right  rotor 

F/A  cyclic  pitch  to  rudder  position  (deg/deg) 

(53)  Coefficient  for  linking  left  rotor 

F/A  cyclic  pitch  to  -udder  position  (deg/deg) 

(54)  Coefficient  for  linking  forward 

rotor  lateral  cyclic  pitch  to  rudder 
position  (deg/deg) 

(55)  Coefficient  for  linking  aft  rotor  (deg/deg) 
lateral  cyclic  pitch  to  rudder  position 

(56) 


(in) 

(deg) 

(deg) 

(deg) 


CARD  99 


(57)  j  Coefficients  for  tail  rotor  col-  (deg) 

(58)  |  lective  as  a  function  of  rudder 

(59)  )  angle  if  XCON(48)/0  (/deg) 

(60)  Coefficient  for  modifying  XCON(36)  and 
XCON(37)  with  fa 

(61)  )  Coefficients  for  modifying  XC0N(50) 

(62)  j  and  X00N(51)  with /8m  (deg) 

(63)  Mast:  tilt  -  elevator  incidence  coupling 

factor  (deg/deg) 
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12 .  Flight  Con s tan ts  Group 


CARD  101 
XFC 


(1)  Forward  velocity 

(2)  Lateral  velocity 

(3)  Rate  of  climb 

(4)  Altitude 

(5)  Euler  angle  yaw 

(6)  Euler  angle  pitch 

(7)  Euler  angle  roll 


CARD  102 

(8)  Collective  stick  position 

(9)  F/A  cyclic  stick  position 

(10)  Lateral  cyclic  stick  position 

(11)  Pedal  position 

(12)  g  level 

(13)  Bank  angle 

(14)  Turn  radius 


CARD  103  i 

(15)  Main  rotor  F/A  flapping  angle 

(16)  Main  rotor  lateral  flapping  angle 

(17)  Tail  rotor  F/A  flapping  angle 

(18)  Tail  rotor  lateral  flapping  angle 

(19)  Main  rotor  thrust 

(20)  Tail  rotor  thrust 

(21)  Trim  type  indicator 

CARD  104 

(22)  Initial  approximation  control 

(23)  Fold  indicator 

(24) 

(25)  Engine  RPM 

(26)  Maximum  engine  horsepower 

(27)  Speed  of  sound 

(28)  Sigma-prime 


(kts) 

(kts) 

(ft/sec) 

(ft) 

(deg) 

(deg) 

(deg) 


(%) 

(%) 

(%) 

(%) 

(deg) 

(ft) 


(deg) 

(deg) 

(deg) 

(deg) 

(lbs) 

(lbs) 


(ft/sec) 


13.  Weapons  Group 


CARD  105 

XGN  (1)  Station  line 

(2)  Buttline  Location  of  weapon 

(3)  Waterline 

(4)  Azimuth 

( 5 )  Elevation 

(6) 

(7) 


(in) 

(in) 

(in) 

(deg) 

(deg) 


19 


14.  Allowable  Error  Group 


CARD  111 

XER 

(1) 

Allowable  error 

in 

(2) 

Allowable  error 
balance 

in 

(3) 

Allowable  error 

in 

(4) 

Allowable  error 
moment  balance 

in 

(5) 

Allowable  error 
balance 

in 

(6) 

Allowable  error 
moment  balance 

in 

(7) 

Allowable  error 
moment  balance 

in 

F/A  force  balance  (lbs) 
lateral  force 


(lbs) 

vertical  force  balance(lbs) 
pitching  and  yawing 

(ft-lb) 


rolling  moment 


(ft-lb) 


main  rotor  flapping 


(ft-lb) 

tail  rotor  flapping 


(ft-lb) 
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Iteration  Limits  Grout 


CARD  121 
XI T 


CARD  122 


CARD  123 


(10) 

(11) 

(12) 

(13) 

(14) 


(15) 


Iteration  limit  for  TRIM 

Iteration  limit  for  rotor  flapping  moment 

balance  in  rotor  analysis 

Partial  derivative  use  indicator 

Partial  derivative  increment  for  STAB 

Trim  output  selector 

Induced  velocity  change  limiter 


Minimum  value  for  main  rotor  flapping 

angle  correction  limit  (deg) 

Minimum  limit  for  tail  rotor  flapping 

angle  correction  limit  (deg) 

Maximum  value  for  use  of  variable 

damper  for  main  rotor  (ft-lb) 

Maximum  value  for  use  of  variable 

damper  for  tail  rotor  (ft-lb) 

Starting  value  for  TRIM  correction 

limit  (deg) 

Minimum  value  for  TRIM  correction 

limit  (deg) 

Maximum  value  for  use  of  vari-  (lb  or  ft-lb) 

able  damper  in  TRIM 


Euler  angle  iteration  selector 
for  TRIM 


16.  Stability  Times  Group 

CARDS  131-132 
TSTAB(l) 

*  Maneuver  times  at  which  stability 

analysis  is  to  be  performed. 


TSTABC14) 


Airload  Printout  Timas  Group 

CARDS  133-134 
TAIR(l) 

[  Maneuver  times  at  which  blade 

.  element  data  is  to  be  printed  out. 

TAIR(14) 


(  sec) 
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18.  rime  Card 
CARD  201 

(1)  Start  time  of  maneuver 

(2)  First  time  increment 

(3)  Time  to  stop  using  first  increment 

(4)  Second  time  increment 

(5)  Time  to  stop  using  second  increment 

(6)  Time  to  stop  the  maneuver 

(7) 


(  sec) 
(  sec) 
(  sec) 
(sec) 
(  sec) 
(  sec) 


19.  Variations  from  Stable  Flight 


CARD  211 

Col.  1 
Col.  2-5 
Col.  11  -  20 
Col.  21  -  30 
Col.  31  -  40 
Col.  41  -  50 
Col.  51  -  60 


Indicator  for  last  card  of  this  group 
J,  variation  selector 

Inputs  which  define  the  variations 
for  each  value  of  J  in  6F10.0 
format 


B.  USER1 S  GUIDE  TO  INPUT  FORMAT 


I.  Identification  and  Flow  Control  Group 
CARD  01 

NPART  is  the  primary  variable  on  this  card.  Uses 
and  meanings  of  the  other  variables  on  this  card 
depend  upon  the  value  of  NPART  and  are  explained 
under  each  value  of  NPART.  If  a  variable  is  not 
listed,  it  is  not  used  for  that  value  of  NPART. 


NPART  <  1 
NPART  =  1 

NPART  =  2 


NPRINT 


NVARB  =  0 


NVARB  =  1 


Data  error,  job  terminated. 

Obtain  stabilized  flight  condi¬ 
tion  only.  Card  134  is  last 
data  card.  See  Figure  1. 

Obtain  stabilized  flight  condi¬ 
tion,  perform  indicated  time- 
variant  maneuver.  CARD  201 
(time  card)  and  at  least  one 
more  card  of  the  200  series  (J 
cards)  must  be  included.  See 
Figures  5  and  6. 

Determines  frequency  of  print¬ 
out  of  maneuver  data.  Program 
prints  data  showing  initial 
conditions  for  maneuver 
(maneuver  time  t=0)  and  every 
NPRINTth  time  point  thereafter. 

Wagner  and  Buettiker  functions 
(non-steady  aerodynamic  effects) 
inactive. 

Wagner  and  Buettiker  functions 
active  for  first  time  increment 
only. 


NVARB  =  2  Wagner  and  Buettiker  functions 

active  for  second  time  increment 
only. 


NVARC  =  0 


Airload  printout  for  main  rotor 
only. 


NVARC  =  1 


Airload  printout  for  both 
rotors. 


NVARC  =  2 


Airload  printout  for  second 
rotor  only. 


NPART  =  3 

NPRINT 

N SCALE  =  0 
N SCALE  =  1 

NS CAL E  =  2 

N SCALE  =  3 
N SCALE  =  4 
N SCALE  =  5 

NSC ALE  =  6 

NSCALE  =  7 

NVARA 

NVARB 

NVARC 

AL( I ) ,  1=1,2 

AH( I) ,  1=1,2 

NPART  =  4 
NPART  =  5 


Plot  time  history  of  selected 
variables.  (See  Figure  7.) 

Determines  frequency  of  time 
points  to  be  plotted.  Used 
in  same  manner  as  under  NPART=2. 

No  effect. 

Multiplies  scales  for  1st  vari¬ 
able  by  1000. 

Multiplies  scales  for  2nd  vari¬ 
able  by  1000. 

Multiplies  scales  for  1st  and 
2nd  variable  by  1000. 

Multiplies  scales  for  3rd  vari¬ 
able  by  1000. 

Multiplies  scales  for  1st  and 
3rd  variable  by  1000. 

Multiplies  scales  for  2nd  and 
3rd  variable  by  1000. 

Multiplies  scales  for  1st,  2nd 
and  3rd  variable  by  1000. 

First  variable,  from  Plot 
index,  Table  I. 

Second  variable,  from  Plot 
index,  Table  I . 

Third  variable,  from  Plot 
index,  Table  I . 

,3  Lower  scale  limit  for  variable 
indicated  by  subscript. 

,3  Upper  scale  limit  for  variable 
indicated  by  subscript. 

Not  used  at  present. 

Not  used  at  present. 


NPART  =  6 
NPART  =  7 

NPART  =  8 

NSCALE  =  0 


NSCALE  =  1 


NSCALE  =  2 


NPAP.T  =  10 


NVARA  =  0 
NVARA  =  1 


Not  used  at  present. 

Stability  analysis.  See  Figures 
3  and  4. 

Time  history  tape  file  opera¬ 
tion.  See  Figures  8  and  9. 

Transfer  time  history  from  plot 
disk  to  storage  tape.  This 
path  allows  subsequent  time  his¬ 
tories  to  be  put  on  same  storage 
tape  if  all  the  time  histories 
are  created  under  the  same  job 
card. 

Transfer  time  history  from  plot 
disk  to  storage  tape.  No  more 
time  histories  may  be  stored  on 
this  particular  tape. 

Transfer  time  history  from  file 
tape  to  plot  disk. 

Parameter  sweep  (trim  or  sta¬ 
bility  derivatives  only).  See 
Figures  2  and  4. 

Trim  only. 

Trim  and  stability  derivatives. 

The  card  or  cards  following 
have  the  new  values  of  the 
sweep  parameters  in  NAMELIST 
input.  The  first  column  of  the 
first  card  must  be  blank.  In 
Columns  2-8,  put:  &CHANGE 
column  9  is  blank.  Starting  in 
column  10  and  continuing  on  as 
many  cards  as  necessary  are  the 
parameters  to  be  changed  and 
their  values.  After  the  last 
one,  have  either  a  comma  or  a 
blank  followed  by  &END.  Any 
variable  on  CARDS  11-123  except 
XFC(5)  through  XFC(ll)  and 
XFC(1?5)  through  XFCC20),  may 
be  changed  by  referring  to  it 
by  its  group  symbol  and  word 
number.  The  excepted  data  are 


NPART  = 


NVARA 
AL(  1 ) 
NVARB 

AL(2) 


the  trimmed  values  from  the 
previous  case. 

Example: 

&CHANGE  XFS(1)=9500.,XMR(11)=5., 
XFC(3)=-50. ,XER=5*1000. ,100. , 
50., YEL(  17 ) = . 1 08 ,&END 

The  above  information  will  cause 
a  case  to  be  run  for  which 
XFC(5 )  through  XFC(ll)  and 
XFC(15)  through  XFC(20)  will 
have  the  trim  values  of  the 
preceding  case,  gross  weight 
will  be  changed  to  9500.  lbs., 
main  rotor  mast  tilt  angle 
changed  to  5.  degrees,  rate  of 
climb  changed  to  -50.  ft/sec, 
the  allowable  force  and  moment 
unbalances  on  the  fuselage 
changed  to  1000  lbs.  or  ft- lbs, 
the  allowable  main  rotor  flap¬ 
ping  moment  unbalance  changed  to 
100  ft-lb,  the  allowable  tail 
rotor  flapping  moment  unbalance 
changed  to  50  ft-lb  and  the 
elevator  lift  curve  slope  for 
M^gh^number  *  0  changed  to 

11  Least  squares  curve  fit  of  time 

history  data.  Variables  which 
can  be  fit  and  their  code  num¬ 
bers  are  found  in  Table  I.  See 
Figure  10. 

Number  of  curves  to  be  fit. 

Frequency  in  cps. 

Number  of  curves  to  be  used  as 
reference  curves  in  computing 
amplitude  ratios  and  phase 
angle  differences. 

Number  of  curves  to  be  expressed 
as  linear  combinations  of  two 
other  curves.  (This  number  must 
have  a  decimal  point  punched.) 


29 


NVARC  Number  of  maneuver  time  points 

to  be  skipped  before  curve  fit 
begins . 

All  cards  in  an  NPART=11  data 
deck,  except  the  NPART  card, 
have  a  1415  format. 

Next  Card(s)  (110X  series)  Code  numbers  of 

curves  to  be  fit  (NVARA  numbers) 

Next  Card(s)  (111X  series) 

Column  1-5:  NX 

Column  6-10:  Code  number  of 
reference  variable  for  computing 
amplitude  ratios  and  phase  angle 
differences. 

Column  11-15  NX  code  numbers 

16-20  of  variables  to 

etc.  be  compared  to 

the  reference 
variable. 

There  are  NVARB  sets  of  cards 
of  this  type. 

Next  Card(s)  (112X  series)  These  cards 

select  the  variables  to  be  ex¬ 
pressed  as  a  linear  combination 
of  two  other  selected  variables. 

A=KB*B+KC*C+D 

Column  1-5:  Code  number  for 

variable  A. 

Column  6-10:  Code  number  for 
variable  B. 

Column  11-15:  Code  number  for 
variable  C. 

The  constants  KB,KC,  and  D  are 
then  determined  by  the  program. 

There  are  AL(2)  cards  of  this 
type. 


30 


CARD  02 


Do  not  use  any  blank  cards  in 
the  NPART  =  11  group. 


IPSN 


CARD  02 
CARD  03 
CARD  04 


Problem 

Printed 


series  number  for  identification  purposes 
in  output  headings. 


Identifying 

Comments 


Printed  in  output  headings; 
used  as  problem  identification. 


2.  Aerodynamic  Inputs 

The  aerodynamic  inputs  are  used  by  CLCD  to  compute  the 
coefficients  of  lift,  Cl,  and  drag,  Cj),  as  a  function  of 
angle  of  attack,  a  ,  and  Mach  number,  M.  In  the  following 
discussion,  YXX(I)  refers  to  the  Ith  aerodynamic  input  for  the 
appropriate  aerodynamic  surface.  Except  for  sections  spec¬ 
ifically  labeled  otherwise,  formulas  and  procedures  apply  to 
all  aerodynamic  surfaces,  i.e.,  rotors,  wings,  elevators,  or 
fins. 

If  the  lift  curve  slope  for  zero  Mach  number,  YXX(17),  is 
input  as  a  zero  value,  the  program  will  compute  Cl  and  Cp  from 
data  tables  for  a  64A210  airfoil.  If  YXX(17)  has  a  non-zero 
value,  Cl  and  Cd  will  be  computed  from  the  analytical/empir- 
ical  equations  described  in  Section  III  B,  Volume  I.  The 
airfoil  section  characteristics  included  in  the  aerodynamic 
inputs  are  for  a  two-dimensional  symmetric  section  at  zero 
Mach  number.  These  input  values  are  corrected  by  the  CLCD  and 
YFIX  subroutines  to  obtain  the  necessary  three-dimensional 
characteristics.  The  Cl  and  Cd  values  computed  from  these 
three-dimensional  characteristics  are  for  an  unswept  symmet¬ 
rical  surface. 


CARD  A 


Drag  divergence  Mach  number  locates  the  boundary 
between  the  subsonic  and  transonic  regions.  The 
program  requires  a  value  less  than  1. 

Mach  number  for  lower  boundary  of  supersonic  region 
locates  the  boundary  between  the  transonic  and 
supersonic  regions.  The  program  requires  a  value 
greater  than 


fi  +  (i  -  ^rxxa)2)(7X^r7T) 


CARDS  A-B 


at  the  stall  point  = 


YXX(3)+YXX(4)M+YXX(5)M2+YXX(6)M3  for 

normal  flow 


YXX(  7  )  +YXX(  8  )M+YXX(  9  )M2  +YXX(  10  )M3f  or 


reversed  flow 


CARD  B 


Tail  boom  bending  coefficient  reduces  Cl,  on  the 
fin  and  elevator  only,  by  the  formula 


CL  =  CL/(1  +  YXX(Il)qtCL/a) 


where  qt  is  the  dynamic  pressure  on  the  appropriate 
tail  surface. 

CARDS  B-C 

For  non-divergent  drag, 
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3.  Fuselage  Grou 


CARD  11 


Gross  weight  includes  pylons,  rotors,  fuel,  etc. 

Fuselage  aerodynamic  center  defines  the  point  of 
application  of  body  lift,  drag,  and  side  force  so 
that  moments  produced  by  them  may  be  calculated. 

C.  G.  location  is  for  whole  ship.  If  a  conversion 
or  fold  maneuver  is  run  and  C.G.  shift  with  mast 
tilt  or  fold  is  desired,  these  inputs  should  be  for 
0  degrees  mast  tilt  and  rotors  unfolded.  (See  CARD 
15,  pylon  weight  and  C.G.) 

CARD  13-14  ,  9 

Let  q  =  1/2  PV 

a  =  Fuselage  pitch  angle  of  attack 

fi  =  Fuselage  sideslip 

Pitch  moment=qf ( XFS(15)  +  XFS(16)a) 

Yaw  moment  =  q'(XFS(17)  +  XFS(18)jS) 

Other  contributions  to  fuselage  moments  will  come 
from  fuselage  and  pylon  forces  if  the  forces  are  not 
applied  at  the  center  of  gravity. 

Forces  calculated  below  are  in  wind  reference. 

Lift  =  q'(XFS(20)  +  XFS(21)a) 

Drag  =  q'(XFS(22)  +  XFS(23)a  +  XFS(24)a2 

+  XFS(25)/32)  +  pylon  drag 

Side  Force  =  q’  (XFS(26)+XFS(2  7)/3+XFS(28)/82) 

CARD  15 

Pylon  weight  is  the  total  weight  of  both  pylons, 
hubs,  rotors,  etc.,  which  contribute  to  C.G.  shift 
with  mast  tilt.  If  the  mast  tilt  angle  is  to  remain 
constant  during  the  run  and  the  C.G.  inputs  of  CARD 
11  are  the  locations  of  the  ship  C.G.  for  the  input 
mast  tilts,  then  pylon  weight  should  be  input  as  0. 
Otherwise,  a  shift  from  the  C.G.  inputs  will  be  cal¬ 
culated  as  explained  below. 


Pylon  C.G.  inputs  are  intended  to  locate  the  movable 
weight  (pylons,  etc.)  for  0  degree  mast  tilt.  Zero 
mast  tilt  means  mast  vertical  for  all  but  a  side 
mount  type  tail  rotor.  For  a  tail  rotor,  zero  mast 
tilt  is  directly  to  the  pilot’s  left. 

Pylon  weight  and  C.G.  location  are  necessary  only  in 
conversion  maneuvers.  However,  they  are  used  in  all 
cases  according  to  the  following  formulas. 

The  shift  of  the  C.G.  station,  ASTA,  and  waterline, 
AWL,  due  to  mast  tilt  are  given  by: 

ASTA  =  Z  sin£m  +  X(1  -  cos^) 

AWL  =  Z  (1  -  cos/S^)  -  X  sin/3m 

Where  is  mast  tilt  angle, 

X  =  [XFS(29)/XFS(1)|  [XMR(8)-XFS( 30)] 

Z  =  [xFS( 29)/XFS(l)]  [XMR( 10)-XFS( 31 )] 

It  is  assumed  here  that  the  C.G.  buttline  does  not 
change  with  mast* tilt  and  that  the  mast  tilt  angles 
of  both  rotors  are  equal  when  this  is  used. 
Essentially,  this  means  that  use  of  this  option 
assumes  that  a  side-by-side  configuration  is  being 
simulated. 

Pylon  differential  flat  plate  drag  area  is  defined 
as  total  ship  flat  plate  drag  area  at  0  degrees  mast 
tilt  minus  total  ship  flat  plate  dra^  area  at  90 
degrees  mast  tilt.  The  drag  from  this  area  is 
applied  as  a  function  of  the  cube  of  the  cosine  of 
the  mast  tilt  angle.  The  drag  from  this  area  is 
applied  at  the  pylon  aerodynamic  center  which  is 
assumed  to  be  on  the  mast  at  some  distance  from  the 
mast  pivot  point  XFS(33),  and  thus  causes  a  pitching 
moment  on  the  fuselage. 


4.  Main  Rotor  Group 
CARD  21 

The  biade  chord  is  assumed  to  be  constant  over  the 
radius. 

The  coning  spring  constant  gives  the  amount  of  re¬ 
storing  moment  generated  by  one  degree  of  coning. 

A  reasonable  value  for  the  coning  spring  constant, 
Kc,  may  be  found  from  the  following  equation. 

2  2  2 

Kc  =  -  Or )  -  elT  ^  r^m(r^)Ar^ 

i  =  l 

where 

1^  =  Mass  moment  of  inertia  for  one  blade. 

<o  =  Natural  frequency  of  the  first  collective  mode, 
ft  =  Rotor  rotational  speed. 

e  =  Distance  from  the  shaft  to  the  flapping  hinge, 
r^  =  Radius  of  ith  radial  station. 

m(r^)  =  Mass  per  unit  length  at  ith  radial  station. 
Ar^  =  Radius/20 

Flapping  stop  is  the  maximum  angle  a  blade  can  flap 
down  relative  to  the  mast.  The  normal  input  is  neg¬ 
ative.  An  input  equal  to  the  precone  angle  gives  a 
rigid  rotor. 

CARD  22 

Mast  tilt  angle  is  positive  for  forward  tilt  of  the 
rotor. 

Mast  length  is  the  distance  from  the  shaft  pivot 
point  to  the  hub. 

Fold  rotor  C.G.  shift  switch=0  gives  no  C.G.  shift 
when  the  rotors  are  folded  or  unfolded.  ^0  does 
give  shift  in  C.G. 

CARD  23 

The  virtual  hinge  point,  XMR(15),  is  used  only  in 
calculating  the  mass  moments  of  inertia  of  the 
rotor.  Only  those  rotor  segments  outboard  of  the 
virtual  hinge  point  are  included  in  the  mass  moments 
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of  inertia.  If  the  real  flapping  hinge  offset, 
XMR(16),  is  non- zero,  then  the  virtual  hinge  point 
is  set  equal  to  the  flapping  hinge  offset. 

The  focused  pylon  coefficient  was  included  to  con¬ 
trol  H-force  and  Y-force  for  a  tilting  prop-rotor. 
Changes  are  made  in  the  swashplate  angles  propor¬ 
tional  to  the  inplane  forces.  These  may  be  expres¬ 
sed  in  equation  form  as 

**1 

and  ABL  =  -  ^  Y 

where 

AA.  is  the  change  in  the  fore  and  aft  swashplate 
angle; 

Abl  is  the  change  in  the  lateral  swashplate  angle; 

H  is  the  H-force; 

Y  is  the  Y-force; 

K  is  the  torsional  spring  rate  which  couples  the 
swashplate  to  the  inplane  forces;  and 

L  is  the  moment  arm  through  which  the  inplane 
forces  are  linked  to  the  swashplate. 

The  ratio  L/K  is  the  program  input  XMR(19), 
focused  pylon  coefficient. 

Flapping  spring  rates  are  the  effective  rates,  per 
blade,  for  a  rigid  hub. 

Hub  extent  is  checked  against  the  radius  of  each 
radial  station.  When  the  radius  is  less  than  the 
hub  extent,  Cl=0,  Cj)=Drag  coefficient  for  hub, 
while  area  is  that  of  the  blade.  That  is,  the  blade 
starts  at  the  center  of  the  rotor  (and  center  of  the 
hub)  but  does  not  generate  lift  inboard  of  the  outer¬ 
most  radial  station  which  is  in  the  hub. 


CARD  24 

Drag  coefficient  for  hub  is  discussed  with  hub  ex¬ 
tent  in  CARD  23.  Radial  flow  factor:  <0  gives  no 
rad:'al  flow;  =0  gives  radial  drag  but  not  radial 
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lift;  >0  gives  both  lift  and  drag  radially.  Normal 
input  is  1/YMR(17). 

The  coefficient  for  tip  vortex  effect  modifies  the 
induced  velpcity  distribution  on  the  outboard  30% 
of  the  rotor  blade  to  simulate  the  effect  of  shed 
tip  vortices.  The  simulation  gives  improved  air¬ 
load  calculations  in  the  low  speed  range.  However, 
power  and  other  performance  values  are  not  affected 
significantly.  Rotor  bending  moments  computed  by 
another  version  of  this  program  showed  improved 
correlation  with  test  data  when  a  value  of  ten  (10) 
was  used  for  this  coefficient.  If  the  input  is  zero 
the  effect  is  removed. 

CARD  25  -  27 

Blade  twist  at  each  radial  station  (up  to  20)  from 
tip  to  root,  washout  negative.  The  usual  set  of 
twist  numbers  is  a  sequence  of  negative  numbers 
decreasing  in  magnitude. 

Angle  of  attack  at  station  j  is  computed  by 


'*j  =  0Q  +  +  TWIST  j 


where  station  1  is  the  tip. 

CARD  28  -  2A 

Blade  weight  distribution  is  lbs/in  for  each  of 
segments  from  root  to  tip  with  the  21st  number  t 
tip  weight  in  pounds. 
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5 .  Tail  Rotor  Group 

See  Main  Rotor  Group  for  all  except: 
CARD  32 


CARD  3U 


Tail  rotor  indicator  =  0.  means  single  rotor  heli¬ 
copter;  =  1.  means  tandem  or  side-by-side  ship. 

Taii  rotor  sidewash  coefficient,  K,  is  used  to  simu- 

late  the  effect  of  the  fuselage  on  the  wind  vector 
as  follows: 


VT  =  Vp  Cl.  -  K) 

where  VF  and  V?  are  the  lateral  components  of  the 

W1Hd4.ueC^°^  1T1  body  reference,  felt  by  the  fuselage 
and  the  tail  rotor,  respectively. 


6.  Wing  Group 


CARD  41 

Location  of  center  of  pressure  of  right  wing  is  for 
the  purpose  of  locating  the  point  of  applications 
of  lift  and  drag  to  compute  their  moments.  The 
wings  are  assumed  to  be  symmetric. 

CARD  42 

Main  rotor  induced  velocity  factor  is  that  fraction 
of  main  rotor  induced  velocity  which  hits  the  wing 
in  the  body  vertical  plane. 

The  deflection  of  the  centerline  of  wing  wake  from 
downwind  =  XWG(9)CL. 

n  =  XWG(10)(CD)1/2cos2(ffD/2h)/( f  +  .3) 

where  D  is  the  distance  from  the  trailing  edge  of 
the  wing  to  the  leading  edge  of  the  elevator  in  wing 
chords,  h  is  the  half  width  of  the  wing  wake  f  wing 
chords  down  the  wake  centerline. 

CARD  42  -  43 

Coefficients  for  XWG( 12)  through  XWG(21)  are  deter¬ 
mined  from  Etkin,  Reference  1,  pp.  486-495,  and 
used  in  the  manner  set  forth  below.  These  variables 
affect  lateral  stability  and  trim  characteristics 
for  non-zero  sideslip. 

Define : 

B  =  Wing  span 
S  =  Wing  area 
£  =  Sideslip 

a  =  Wing  angle  of  attack 
p  =  Roll  rate  of  fuselage 
r  =  Yaw  rate  of  fuselage 
L  =  Roll  moment  of  wings 
N  =  Yaw  moment  of  wings 
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Then 


L  =  F[/3(XWG(12)+XWG(13)CL)+t*(XWG(14)rCL+XWG(15)p)] 
and 

N  =  f|/3(XWG(16)+XWG(17)Gl2) 

+  t*|r(XWG(18)CL2+XWG(19)dCD/da) 

+  p(XWG(20)CL+XWG(2L)(YWG(L3) 

+  2aYWG(14)YWG(18)/(3+YWG(18))))|! 


7 .  Elevator  Groui 


CARD  SI 

Location  of  the  center  of  pressure  is  the  point  of 
application  of  lift  and  drag  to  determine  moments. 
XEL(5)  is  positive  leading  edge  up. 

CARD  52 


is  the  induced  velocity  of  the  main  rotor 
V  is  the  airspeed 

(  -  XEL(11)Cl  where  (  is  the  change  in  elevator 

angle  of  attack,  Cl  is  the  average  of  the  Cl 
on  both  wings 


8.  Fin/Rudder  Groui 
CARD  61 


Location  of  the  center  of  pressure  is  the  point  of 
application  of  lift  and  drag  to  determine  moments. 

Geometric  angle  is  considered  positive  if  it  tends 
to  increase  the  lateral  force  which  is  positive  to 
the  pilot’s  right. 

Induced  velocity  factor  is  for  the  portion  of  the 
tail  rotor  induced  velocity  which  affects  the  fin. 


Sidewash  coefficient  is  used  to  simulate  the 
effect  of  the  fuselage  on  the  wind  vector 


Vf  =  VB(1-  ‘ 


where  Vf  and  Vb  are  the  lateral  components  of  the 
wind  vector,  in  body  reference,  felt  by  the  fin  and 
fuselage, respectively. 
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9.  Jet  Group 
CARD  71 


Number  of  jets  should  be  0.,  1.,  or  2.  If  number  of 
jets  =  1.,  left  jet  thrust  is  not  changed  from  the 
input  value  by  either  movement  of  the  collective 
stick  (if  XC0N( 12)^0)  or  by  maneuver  inputs. 

CARD  71  -  72 

Location  of  right  jet  is  the  point  of  application  of 
thrust  with  the  angles  defining  the  direction  of 
thrust  with  respect  to  the  fuselage. 

A  positive  yaw  angle  gives  the  right  jet  a  force 
component  in  body  reference  which  is  positive  to 
the  pilot’s  right.  A  positive  pitch  angle  gives 
the  right  jet  a  force  component  in  body  reference 
which  is  positive  up.  The  left  jet  is  assumed  to  be 
symmetric  in  position  and  attitude  with  the  right 
jet. 
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10.  Bobweight  Group 
CARD  81  -  83 

For  no  bobweight,  set  rj  =  0. 

For  a  bobweight,  set  rj  >  0. 

The  parameters  listed  in  the  Input  Format  are  for  a 
special  bobweight  design  and  are  defined  in  Refer¬ 
ence  2.  A  value  of  »?  =  4  is  assumed  in  that 
analysis. 

The  equation  for  the  basic  case  is 

355  5  +  CXS  +  KjjS  =  -5^3  355  '  max  j 

V  Cg-gp) 

where 

5  =  change  in  collective  pitch  due  to  the  bobweight 
(radians) 


g  =  vertical  acceleration  of  the  fuselage 

The  remaining  symbols  are  described  in  the  input 
list  below: 


CARD  81 
XBW 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 


Effectivity  coefficient 
Spring  constant 
Damping  coefficient 


(deg/in) 
( in- lb/rad) 
(lb-in- sec) 


0. 

0. 

0. 

g 

a 


? 


Preload  (bobweight  not  effective 
g<gp) 


(g’s) 


CARD  82 

(8)  0.  2 
(9)  I  Weight  moment  of  inertia  of  bob-  (lb-in  ) 
weight 
(10)  0. 

(11)  0. 

(12) 

(13) 

(14) 


,•*  . 
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CARD  83 


(15)  1 

(16)  1 

(17)  1 

(18)  1 

(19)  1 

(20)  1 

(21)  1 


11.  Controls  Group 

a.  Collective 
CARD  91 

Lower  limit  =  XCON(  2)+XC0N( 3)£m+XC0N( 4)£m2 
Range  =  XC0N(5)+  [xCON(6)-XCON(5)]/3m/90 

CARD  92 

A  collective  stick  movement  such  that  the  collective 
pitch  of  rotor  1  is  increased  by  one  degree  results 
in  a  collective  change  on  rotor  2  of  XCON(9)  degrees. 

A  stick  movement  of  one  inch  results  in  a  change  of 
wing  angle  of  attack  of  XCON(IO)  degrees,  a  change 
of  elevator  angle  of  attack  of  XCON(ll)  degrees, 
and  a  change  in  the  thrust  of  each  jet  of  XCON(12) 
pounds. 

Maximum  rate  for  PCG  is  the  maximum  rate  which  the 
PCG,  when  active*  can  change  the  root  collective 
pitch. 

RPM  dead  band  is  ARPM/RPMref  where  ARPM  is  the  maxi¬ 
mum  change  in  RPM  which  does  not  activate  the  PCG 
and  RPMref  is  the  input  RPM. 

b.  F/A  Cyclic 
CARD  93 

A  movement  of  the  stick  through  100%  of  its  range 
is  equivalent  to  movement  of  the  stick  through 
XC0N(15)  inches.  During  a  stick  movement  from 
full  aft  to  full  forward,  the  swashplate  F/A  set¬ 
ting  will  move  from  XC0N(16)  degrees  to  X00N(16)  + 
XCON(17)  degrees. 

CARD  94 

A  F/A  cyclic  stick  movement  such  that  the  F/A  cyclic 
pitch  of  the  forward/right  rotor  is  increased  by 
one  degree  results  in  a  collective  pitch  change  on 
the  forward  rotor  of  XCON(22)  degrees,  a  collective 
pitch  change  on  the  aft  rotor  of  XCON(23)  degrees, 
a  F/A  cyclic  pitch  change  on  the  aft/left  rotor 
of  XCON(24)  degrees,  a  lateral  cyclic  pitch  change 
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on  the  right  rotor  of  XCON(25)  degrees,  and  a  lateral 
cyclic  pitch  change  on  the  left  rotor  of  XCON(26) 
degrees.  However,  XCON(25)  and  XCON(26)  are  used 
only  when  XCON(33)->2  and  the  result  of  their  in¬ 
fluence  is  added  to  XCON(32). 

Aa£  =  XCON(27)K  +  XCON(28)K2  where  is  the  change 
in  geometric  angle  of  incidence  of  the  elevator  and 
K  is  the  displacement,  in  inches,  of  the  stick  from 
neutral . 

c.  Lateral  Cyclic 

CARD  95 

XCON( 29 ) ,  XCON(30),  and  XC0N(31)  are  the  lateral 
cyclic  equivalents  of  XC0N(15),  XC0N(16),  XC0N(17). 

CARD  96 

A  lateral  cyclic  stick  movement  such  that  the 
lateral  cyclic  pitch  of  the  forward/right  rotor  is 
increased  by  one  degree  results  in  a  collective 
pitch  change  on  the  right  rotor  of  XC0N(36)  degrees 
a  collective  pitch  change  on  the  left  rotor  of  XCON 
(37)  degrees,  a  F/A  cyclic  pitch  change  on  the  right 
rotor  of  XCON(38)  degrees,  a  F/A  cyclic  pitch  change 
on  the  left  rotor  of  X0CN(39)  degrees,  and  a  lateral 
cyclic  pitch  change  on  the  aft/left  rotor  of  XCON(40) 
degrees. 

A  lateral  cyclic  stick  movement  of  one  inch  to  the 
right  results  in  an  increase  in  the  left  wing  angle 
of  attack  of  XC0N(41)  degrees  and  an  equal  decrease 
in  the  right  wing  angle  of  attack. 

d.  Pedal  and  Mast  Tilt 


CARD  97 


X00N(43),  XC0N(44),  XCON(45)  are  the  pedal  equiva¬ 
lents  of  XC0N(15),  XC0N(16),  XC0N(17)  and  apply  to 
tail  rotor  collective  if  XC0N(48)=0,  or  apply  to 
the  rudder  if  XCON(48)/0. 

For  a  tandem  or  side-by-side  configuration  which 
does  not  have  either  a  tail  rotor  or  a  rudder,  it 
is  necessary  to  input  numbers  for  a  fictitious 
rudder  (i.e. ,  XC0N(48)=0)  in  XC0N(43),  XCON(44), 
and  XC0N(45)  in  order  to  define  the  inputs  on 
CARD  98.  It  is  best  to  construct  this  fictitious 
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rudder  so  that  at  least  one  of  the  inputs  on  CARD 
98  will  be  equal  to  1. 

Rudder  lock  indicator  =0.  means  fin  is  fixed  and 
does  not  move. 

Rudder  lock  indicator  =  1.  means  the  fin  angle  of 
attack  is  directly  proportional  to  pedal  deviation 
from  neutral. 


CARD  98 


A  pedal  movement,  to  the  left,  such  that  the  rudder 
angle  or  tail  rotor  collective  is  increased  by  one 
degree  results  in  a  collective  pitch  change  on  the 
right  rotor  of  XCON(50)  degrees,  a  collective  pitch 
change  on  the  left  rotor  of  XC0N(51)  degrees,  a  F/A 
cyclic  pitch  change  on  the  right  rotor  of  XCON(52) 
degrees,  a  F/A  cyclic  pitch  change  on  the  left  rotor 
of  XCON(53)  degrees,  a  lateral  cyclic  pitch  change 
on  the  forward  rotor  of  X00N(54)  degrees,  and  a 
lateral  cyclic  pitch  change  on  the  aft  rotor  of 
X00N(55)  degrees. 


CARD  99 


e  =  XCON(  57  )  +XC0N(  58)P+X00N(  59  )P3  where  6  is  the 
collective  pitch  on  the  tail  rotor  of  a  single  rotor 
helicopter  and  P  is  the  deviation,  in  inches,  of 
the  pedal  from  neutral. 

XC0N(36)  and  XCON(37)  are  modified  by  cos(XC0N(60) 

Am)  • 

XCON(50  and  XC0N(51)  are  modified  by  sir(XC0N(61) 

+  XCONC62)). 

An  increase  in  mast  tilt  angle,  0m,  of  one  degree 
resul ts  in  an  increase  in  elevator  angle  of  attack 
of  XCON(63)  degrees. 

The  following  inputs  are  modified  by  cos  An: 
XCONC22),  XCdN( 23) ,  XC0N(24),  XCON(25),  XC0N(26) , 
XCON(40),  XCONC52),  and  XC0N(53)f  while  XC0N(38) 
and  XC0N(39)  are  modified  by  sin  0m. 

Inputs  in  the  Controls  Group  are  used  only  where 
appropriate  as  defined  by  the  configuration  control 
linkage  tables,  Table  11  and  Table  111. 
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If  the  following  control  inputs  are  zero,  they  are 
reset  to  the  values  indicated. 


XCON(L)  -  100. 

XC0N( 5 )  -  100. 

XC0N( 6)  -  100. 

XCON ( 9 )  -  I. 

XC0N( 15)  -  100. 

XG0NCI6)  -  -50. 

XC0NCI7)  -  100. 

XC0N( 24)  -  I. 

XG0N( 29 )  -  100. 

XCON( 30 )  -  -50. 

XG0N(3I)  -  100. 

XC0N( 40)  -  I. 

XCONC43)  -  100. 

XG0N( 44)  -  -50. 

XC0N( 45 )  -  100. 

Control  deviation,  in  inches,  from  neutral  is 
computed  by 

K  =  (K^/IOO  +  K3)  K4/K2 

where 

is  the  control  position  in  percent 

K2  is  the  range  of  the  control  surface  in  degrees 

K~  is  the  lower  limit  of  the  control  surface  in 
degrees 

is  the  range  of  the  control  in  inches 


CARD  101 


CARD  102 


CARDS  102 


CARD  104 


The  input  velocities  are  with  respect  to  the  ground. 

The  Euler  Angles  are  the  angles  between  the  coordi¬ 
nate  system  attached  to  the  body  and  the  ground 
coordinate  system.  Yaw  is  positive  nose  right; 
pitch  is  positive  nose  up;  roll  is  positive  down 
right. 

Altitude  is  the  height  above  ground.  It  is  used  in 
the  calculations  for  ground  effect.  If  altitude  is 
negative,  the  program  stops. 


The  program  checks  g-level,  bank  angle,  and  turn 
radius,  in  that  order,  for  a  non-zero  input.  The 
first  non-zero  input  completely  specifies  the  turn 
and  the  program  ignores  the  other  inputs.  During 
the  TRIM  process,  instead  of  seeking  zero  force  and 
moment  values,  the  program  seeks  those  force  and 
moment  values  which,  during  an  undisturbed  maneuver, 
will  cause  the  ro*torcraft  to  fly  in  the  circle 
determined  by  the  input.  A  g-level  input  will  pro¬ 
duce  a  right  turn.  A  left  turn  can  be  produced  by 
a  negative  bank  angle  or  a  negative  turn  radius. 

-  103 

Flapping  angles  on  both  rotors,  control  positions, 
and  Euler  angles  of  pitch  and  roll  are  changed  in 
the  TRIM  subroutine  as  necessary  in  order  to  achieve 
stable  flight  conditions. 

Trim  type  indicator:  =  0.  is  for  either  a  coordinated, 
banked  turn,  as  explained  above  for  CARD  102,  or  a 
Trim  with  no  turn.  =1.  is  for  a  pull-up  or  push¬ 
over  at  the  g-level  input  in  XFC(12). 


Initial  approximation  control  =  0.  means  the  pro¬ 
gram  is  to  make  an  initial  approximation  to  the 
flapping  angles  and  thrust  on  both  rotors,  ignoring 
the  inputs. 

Initial  approximation  control  =  1.  means  the  pro¬ 
gram  is  to  use  the  inputs  for  flapping  on  both  rotors 
and  make  an  initial  approximation  for  thrust  on  both 
rotors,  ignoring  the  thrust  inputs. 
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.✓*  • 

/ 

Initial  approximation  control  =  2.  means  the  pro¬ 
gram  uses  the  inputs  for  the  flapping  angles  and 
thrust  on  both  rotors  as  the  initial  approximations. 

Fold  indicator  =0.  if  the  rotors  are  unfolded  and 
at  the  RFM  specified,*  /0  if  the  rotors  are  stopped 
and  folded.  ^0  is  used  to  trim  for  an  unfold  and 
start  maneuver.  The  data  should  be  set  up  as  if 
the  rotors  were  unfolded  and  at  normal  RFM  except 
for: 

1.  XFC(  23)7*0 

2.  Controls  are  looked  by 

XC0N( 7 ) , XCCN( 18 ) , XC0N( 32 )  ,  XCON( 46 )/0 

3.  Maneuver  input  cards  for 

J=18  and  J=27  have  a  start  time  of  0,  i.e., 
for  J  =  18,  Q0>O 

Sigma-prime  is  the  ratio  of  air  density  to  standard 
day  air  density.  P  =  afP0  where  PD  =  .002378  slugs/ 
ft^.  p  is  then  used  in  the  calculation  of  dynamic 
pressure. 


13.  Weapons  Group 
CARD  105 

Stationline,  buttline  and  waterline  are  used  to 
locate  the  point  of  application  of  the  recoil  force 
of  the  weapon. 

Azimuth  and  elevation  define  the  orientation  of  the 
weapon  with  respect  to  the  fuselage. 

If  both  angles  are  input  positive,  the  gun  is  point¬ 
ing  up  and  to  the  pilot’s  right. 
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15 .  Iteration  Limits  Group 
CARD  121 


If  the  program  has  not  found  a  stable  flight  condi¬ 
tion  after  XIT(l)  iterations,  it  assumes  there  is 
no  stable  flight  condition  and  terminates  the  pro¬ 
blem. 

* 

If  the  iteration  limit  for  the  rotor  flapping  moment 
balance  in  the  rotor  analysis  is  input  as  a  non-zero 
number,  the  TRIM  equations  are  decoupled  into  one 
set  of  six  equations  in  six  unknowns  and  two  sets 
of  two  equations  in  two  unknowns.  This  should  be 
done  only  on  cases  for  which  the  coupled  10  x  10 
system  has  defied  solution  because  run  time  is  in¬ 
creased  by  a  factor  of  3  to  5. 

Partial  derivative  use  indicator:  =  0.  means 
recompute  partial  derivative  matrix  for  each  itera¬ 
tion  in  TRIM.  /  0.  means  compute  the  partial  deri¬ 
vative  matrix  on  iterations  1,  6,  11,  etc.,  in  TRIM. 
During  the  other  iterations,  use  the  last  computed 
matrix. 

Partial  derivative  increment  for  STAB,  XIT(4),  is 
used  in  computing  rate  derivatives.  The  angular 
rate  increment  is  1/10  the  input,  in  radians/sec, 
and  the  linear  rate  increment  is  10  times  the  input, 
in  ft/sec. 

Trim  output  selector,  XIT(5),  =  0.  prints  the  last 
set  of  partial  derivatives  and  forces  and  moments 
in  TRIM  and  a  set  of  variables  for  each  time  point 
in  MANU.  =  1.  prints  each  set  of  partial  derivatives 
and  forces  and  moments  in  TRIM  and  a  set  of  variables 
for  each  time  point  in  MANU. 

The  induced  velocity  change  limiter  is  twice  the 
maximum  amount  the  induced  velocity  is  allowed  to 
change  between  time  points  during  a  maneuver. 


CARD  122 


XIT( 12 ,13,14)  are  the  inputs  for  a  variable  damping 
procedure  in  TRIM.  At  each  iteration,  a  delta  is 
computed  for  each  of  the  10  trim  variables.  These 
deltas  are  then  compared  to  a  correction  limit.  If 
any  of  the  deltas  are  larger  than  the  limit,  they 
are  ratioed  down  such  that  the  largest  one  is  equal 
to  the  limit.  The  starting  value  of  the  limit  is 
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input  as  XIT(12).  If  at  the  start  of  any  iteration 
the  force  and  moment  errors  are  all  less  in  absolute 
value  than  XIT(14),  then  the  limit  is  halved  with 
the  restriction  that  the  limit  is  never  less  than 
XIT(  13) .  If  XIT(12)<0.5  deg  or  >10  deg,  the  pro¬ 
gram  resets  it  to  L  degree.  If  XIT(13)<.05  deg  or 
>1.0  deg,  the  program  resets  it  to  .05  deg.  If 
XIT(14)<40  XER(I),  the  program  resets  it  to  40  XER(l). 


The  partial  derivative  increment  for  TRIM  is  com¬ 
puted  from  the  correction  limit.  XIT(8,9,10,II) 
define  variable  dampers  for  balancing  the  rotor 
flapping  moments  in  the  rotor  analysis.  The  starting 
value  for  the  correction  limit  is  computed  from  the 
minimum  value  input. 


CARD  123 


The  Euler  angle  iteration  selector  for  TRIM=1  for 
iterating  on  pi*-ch  and  roll;  -2  for  iterating  on 
pitch  and  yaw;  .  l  or  2  gives  =  1. 
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16.  Time  Card 


CARD  201 

If  a  secondary  time  increment  is  not  desired,  only 
start  time,  first  time  increment,  and  time  to  stop 
using  first  time  increment  are  necessary  inputs. 

If  a  secondary  time  increment  is  desired  but  a 
second  use  of  the  first  time  increment  is  not  de¬ 
sired,  the  time  to  stop  the  maneuver  should  be 
either  blank  or  equal  to  the  time  to  stop  using  the 
second  time  increment. 

This  card  and  subsequent  cards  are  to  be  put  into 
the  data  deck  only  when  running  a  maneuver,  i.e., 
NPART=2  on  CARD  01.  The  time  increment  is  the  At 
in  the  Runge-Kutta  solution  of  the  differential 
equations  of  motion.  It  may  be  desirable  to  chnage 
the  value  oi  At  because  of  changes  in  rotor  speed. 
Provisions  have  been  made  for  such  a  change.  See 
CARD  211,  J=31  for  an  example  of  another  use. 

Unless  specified  otherwise,  all  times  are  relative 
to  the  scale  established  by  the  inputs  on  this  card. 
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17.  Variations  from  Stable  Flight 


CARD  211 


NEXT  =  0  This  is  the  last  card  of  the  211  type. 

NEXT  £  0  Another  card  of  the  211  type  follows. 

J  Type  of  variation,  explained  in  list 

below. 

It  is  pof  }ible  to  have  as  many  as  20  cards  similar 
to  CARD  211.  All  have  the  same  format.  It  is  not 
necessary  to  have  the  J  values  in  order,  in  fact, 
there  may  be  several  cards  with  the  same  value  of  J. 
It  is  necessary  that  NEXT  /  0  on  all  of  these  cards 
except  the  last  one  which  must  have  NEXT  =  0. 

When  the  variations  to  be  made  during  the  maneuver 
have  been  decided,  the  list  below  may  be  used  to 
find  the  proper  values  of  J. 


r  ■% 


j  =  i 


J  =  2 


J  =  3 


J  =  4 
J  =  5 


J  =  6 


J  =  7 
J  =  8 


J  =  9 


i 


is  for  movement  of  collective  stick 

is  for  movement  of  F/A  cyclic  stick 

is  for  movement  of  lateral  cyclic  stick 

is  for  movement  of  pedal 

is  for  the  prop-rotor  collective 
governor 

is  for  folding  rotors  aft  after  tilting 
forward  and  stopping 

is  for  a  change  in  rpm  governor  setting 

is  for  converting  the  tail  rotor  into  a 
pusher  prop.,  i.e.,  change  the  tail 
rotor  mast  tilt  angle  in  the  horizontal 
plane. 

is  for  a  vertical  ramp  gust,  ramp  length 
may  be  0. 


J  =  10 
J  =  11 

J  =  12 
J  =  13 

J  =  14 

J  =  15 
J  =  16 
J  =  17 

J  =  18 
J  =  19 
J  =  20 

J  =  21 

J  =  22 

J  =  23 
J  =  24 
J  =  25 
J  =  26 
J  =  27 

J  =  28 
J  =  29 

J  =  30 
J  =  31 


is  for  a  vertical .sine-squared  gust 

a  horiz°ntal  ramP  gust,  ramp 
length  may  be  0. 

is  for  a  horizontal  sine-squared  gust 

is  for  a  change  in  engine  torque 
supplied  M 

plied*7  3  chatl®e  au*iliary  thrust  su] 

is  for  activation  of  a  "yaw  pilot", 
is  for  weapon  fire 

is  for  change  of  mast  tilt  angle  on 
both  rotors  and  rpm  change 

is  for  rotor  brake 

is  for  pitch  A.S.E. 

V 

is  for  sinusoidal  movement  of  controls 
or  mast 

is  for  a  flat  tracker  for  the  main 
rotor 

rotor*7  3  flat  tracker  for  the  tail 

is  for  rpm  dependent  hubsprings 
not  used 
not  used 
not  used 

stop0r  folding  ™tor.  horizontally  after 

is  for  rpm  dependent  flapping  stops 

is  for  connecting  and  disconnecting 
helicopter  controls  6 

is  for  rotor  moment  balancing  mechanism 

is  for  changing  NPRINT  on  CARD  01 
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After  the  values  of  J  have  been  determined,  the  list 
below  gives  details  for  determining  the  rest  of  the 
inputs . 

J  -  1,2 ,3, 4  (Control  Movements) 


Col.  11-20  Start  time  (sec) 

21-30  Input  rate  1  (in/sec) 

31-40  Stop  time  (sec) 

41-50  Start  time  (sec) 

51-60  Input  rate  2  (in/sec) 

61-70  Stop  time  (sec) 


J  =  5  (Proprotor  Collective  Governor) 

Same  as  for  J  =1,  except  rates  are  in  hp/sec. 
J  =  6  (Folding  Rotors  Aft) 


Col.  11-20  Start  time  (after  12=0)  (sec) 

21-30  Rate  (Positive  to  fold  aft)  (deg/sec) 

31-40  Stop  time  (after  0=0)  (sec) 

41-50  Start  time  (after  Q=0)  (sec) 

51-60  Rate  (Positive  to  fold  aft)  (deg/sec) 

61-70  -Stop  time  (after  S2»0)  (sec) 


J  =  7  (Governor  Rpm  Setting) 

Same  format  as  for  J  =  1,  except  rates  are  in 
rpm/sec- 

J  =  8  (Tail  Rotor  to  Pusher  Prop) 

Same  format  as  for  J  =  1,  except  rates  are  in 
deg/sec. 

J  =  9  (Vertical  Ramp  Gust) 

See  J  =  11 

J  =  11  (Horizontal  Ramp  Gust) 


Col.  11-20  (1)  Starting  distance  (ft) 

21-30  (2)  1st  max.  velocity  positive  (ft/sec) 

down  or  north 

31-40  (3)  1st  ramp  length  (ft) 

41-50  (4)  Distance  gust  is  steady  (ft) 

51-60  (5)  2nd  ramp  length  (ft) 


61-70  (6)  2nd  max.  velocity  (measured  (ft/sec) 

from  first  max  velocity) 
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Ramp  Gust 

J  =  10  (Vertical  Sine-Squared  Gust) 
See  J  =  12  ~~ 

— (Horizontal  Sine-Squared  Gust) 


GoL*  starting  distance 

^1-30  (2)  1st  max  value  (positive 

down  or  north) 

31-40  (3)  1st  gust  length 

Distance  between  gusts 
2nd  Sust  length 
oi-70  (6)  2nd  max  value 


(ft) 

(ft/sec) 

(ft) 
(ft) 
(ft) 
(f t/sec) 


Sine -Squared  Gust 


.  4  is 


J  =  13  (Main  Engine  Torque) 


Col.  11-20 

21-30 

31-40 

41-50 

51-60 


Start  time  for  rotor  torque  (sec) 

supplied  variation 

Ratio  of  torque  desired  to  (f t-lb/f t-lb) 
torque  required  at  trim  point 
Start  time  for  rotor  torque  (sec) 

supplied  recovery  to  torque  required 
( inactive) 

Engine  accel.  lag,  zero  to  full  (sec) 
power 


J  =  14  (Auxiliary  Jet  Thrust) 


Col.  11-20 
21-30 


31-40 

41-50 

51-60 

61-70 


Start  time  for  auxiliary  thrust  (sec) 

variation 

=0  or  1  see  below: 

=2  for  setting  col.  51-60  to  trim 
thrust  and  resetting  col.  21-30  to  1. 
Rate  (lb/sec) 

Stop  time  if  col.  21-30  =  0.  (sec) 

Final  value  of  thrust  if  col.  (lb) 

21-30  =  1. 

=1  for  left  jet;  =2  for  right  jet 


J  =  15  (Yaw  Reactions  by  Pilot) 


Col.  11-20 
21-30 

31-40 

41-50 

51-60 

61-70 


Start  time  (sec) 

Acceptable  sideslip,  not  (deg) 

requiring  correction 

T/R  collective  rate  per  (deg/sec/deg) 
sideslip,  X4 

T/R  collective  rate  per( deg/sec/deg/sec) 

sideslip  rate,  X5 

Ratio  of  acceptable  yaw  velocity 

to  maximum  yaw  velocity 

Stop  time 


If  T/R  collective  is  to  be  changed, 
rate  =  /3X3  +  /3X4 

J  =  16  (Machine  Gun  Fire,  Ramp  Only) 


11-20 

Start  time 

(sec) 

21-30 

Stop  time 

(  sec) 

31-40 

Max  force 

(lb) 

41-50 

Ramp  length 

(  sec) 

51-60 

(inactive) 

61-70 

(inactive) 
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J  =  17  (Mast  Tilt  on  Both  Rotors) 


Col.  11-20 
21-30 
31-40 
41-50 
51-60 

61-70 


Start  time  for  mast  tilt  (sec) 

Rate  of  mast  tilt  (deg/sec) 

Stop  time  for  mast  tilt  (sec) 

( inactive) 

a,  mast  tilt  angle  at  which  rpm  (deg) 
change  is  activated 

ftn  _  ^At  change  in  rpm  in  conver-  (rpm) 
titig  from  airplane  mode  to  heli¬ 
copter  mode 


[fiA  +  (flR 


ft^)  cos 


An  -  “ 


.56  a 


ft  =< 


«H  if  An«  “ 


if  A»>  0 

m 


where:  )8m  is  mast  tilt  angle 

ft  is  current  rotor  rpm 

ftH  is  rotor  rpm  in  helicopter  mode  (0m  =  0°) 
ftA  is  rotor  rpm  in  airplane  mode  (/Sm  =  90°) 
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J  =  18  (Rotor  Brake) 


Col.  11-20  Maximum  brake  torque  (ft-lb) 

21-30  RPM  at  which  brake  engages,  flu  (rpm) 

31-40  Target  azimuth  position  for  stop  (deg) 

41-50  Time  to  stop  applying  brake  (sec) 

J  =  19  (Pitch  A.S.E.) 

Col.  11-20  Start  time  (activate  A.S.E.)  (sec) 

21-30  Cyclic  displacement/pitch  (deg/deg) 

displacement  from  reference 
pitch  angle,  X2 

31-40  Cyclic  displacement/pitch  (sec) 

velocity,  X3 

41-50  Cyclic  displacement/pitch  2 

acceleration,  X4  (sec  ) 

51-60  Reference  pitch  angle,  X5  (deg) 

61-70  Lag  factor  for  mechanism,  X6  (sec) 


Let  5=  change  in  cyclic  pitch  for  one  time 
increment,  At 

value  of  £  for  the  previous  time  increment 
a  =  fuselage  Euler  angle  pitch 
q=  fuselage  pitch  rate 


5  = 


J  =  20  (Sinusoidal  Movement  of  Controls  or  Mast) 


11-20 

Start  time 

(sec) 

21-30 

Frequency 

(cps) 

31-40 

Amplitude  (inches  for  controls, 

deg.  for  mast) 

41-50 

Stop  time 

(sec) 

51-60 

Control  to  be  moved 

1.  is  for  collective  stick 

2.  is  for  F/A  cyclic  stick 

3.  is  for  lateral  cyclic  stick 

4.  is  for  pedal 

5.  is  for  mast 

J  =  21  (Flat  Tracker  for  Main  Rotor) 


See  J  =  22 
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J  =  22  (Flat  Tracker  for  Tail  Rotor) 


Col.  11-20 
21-30 
31-40 
41-50 
51-60 
61-70 


Start  time 

Acceptable  deviation 
Rate 

F/A  reference  angle 
Lateral  reference  angle 


Data  selector  0.  -  use  input  values 
for  reference  angles;  1.  =  use  TRIM 
values  of  flapping  angles  wrt  mast 
as  reference  angles 


=  23  (RPM  Dependent  Hubsprings) 


(  sec) 
(deg) 
(deg/sec) 
(deg) 
(deg) 


Col.  11-20  Rotor  number  (1  or  2) 

21-30  Kg  hubspring  value  in  lower  (ft-lb/deg) 
rpm  range 

31-40  top  of  lower  rpm  range  (rpm) 

41-50  ^2  bottom  of  upper  rpm  range  (rpm) 

Let  Q  be  the  rpm  of  rotor  lt  Kj  be  one  of  the 
following  inputs,  as  appropriate;  XMR(17),  XMR(18), 
XTR(17)  or  XTR(18),  and  Kh  be  the  rpm  dependent  value 
of  the  appropriate  hubspring.  Then 


11-20 

21-30 


Kh  = 


if  a  >a. 


(n^)M 


+  K, 


if 


if  ftL<  Q<ft2 


In  other  words,  the  extreme  values  for  the  hub- 
springs  are  Kb,  and  the  inputs  in  the  rotor  groups 
and  linear  interpolation  is  used  in  the  transition 
region. 

J  =  27  (Horizontal  Fold,  for  Main  Rotor  Only) 

Col.  11-20  Start  time  (sec  after  ft=0) 

21-30  Rate  (deg/sec) 

31-40  Stop  time  (sec  after  ft=0) 

41-50  Blade  number  (each  blade  moves 
independently) 
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J  =  28  (Rpm  Dependent  Flapping  Stops) 

Same  as  for  J=23  except  that  mechanism  affected  is 
flapping  stops  and  Kg  is  in  degrees. 

J  =  29  (Control  Changer  -  to  Lock  or  Unlock 
Swashplate) 


Col.  11-20 

Start  time 

(  sec) 

21-30 

Stop  time 

(  sec) 

31-40 

Indicator:  =0 

if 

start  time  is 

in 

maneuver  sec., 
sec.  after  fl=0 

/o 

if  start  time 

is 

in 

1 

V 

41-50 

Indicator:  =0 

if 

stop  time  is 

in 

* 

maneuver  sec. , 
sec.  after  fi=0 

/o 

if  stop  time 

is 

in 

Note:  If  this  mechanism  is  off  during  a  maneuver, 

swashplate  settings  will  immediately  assume 
the  value  dictated  by  the  control  positions. 
Care  should  be  taken  to  set  the  controls  so 
that  there  are  no  discontinuities. 

51-60  Indicates  which  controls  to  lock  or 

unlock:  1.  is  for  rotor  1  collective; 

2.  is  for  F/A  cyclic;  4.  is  for  lat. 
cyclic;  8.  is  for  rotor  2  collective; 
for  any  combination,  add  the  indicators. 
0.  is  equivalent  to  15^  which  affects 
all  controls. 


J  =  30  (Mechanism  for  Balancing  Main  Rotor  Force 
and  foments  During  riorizontal  Fold) 


Col.  11-20 
21-30 
31-40 
41-50 

51-60 

61-70 


Start  time 
Stop  time 
dZ-force/d  collective 
dF/A  flapping  moment/d 
F/A  cyclic 

d  Lat  flapping  moment/d F/A 
cyclic 

Maximum  rate  of  change  of 


(sec  after  Q=0) 
(sec  after  ft=0) 
(lbs /deg) 
(f  t-lb/deg) 


(ft-lb/deg) 


controls  (collective  and  cyclic) 


(deg/sec) 
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J  =  31  (Changing  Printout  Frequency) 


Col.  11-20 
21-30 
31-40 
41-50 
51-60 
61-70 


Time  to  change  NPRINT 
New  NPRINT 

Time  to  change  NPRINT 
New  NPRINT 

Time  to  change  NPRINT 
New  NPRINT 


(  sec) 
(  sec) 
(  sec) 


NPRINT  must  be  input  as  a  floating  number,  i.e., 
punch  a  decimal  point  on  the  data  card.  The  use  of 
NPRINT  is  as  described  for  CARD  01,  NPART=2 . 

As  an  example  of  the  use  of  this  value  of  J,  as  well 
as  an  example  of  the  use  of  the  provision  for  dif¬ 
ferent  time  increments  on  CARD  201,  consider  the 
following  hypothetical  situation. 

A  maneuver  was  run  in  which  a  pitch  divergence 
occurred.  Analysis  of  the  output  indicated  that 
the  divergence  started  between  3.5  and  3.75  seconds. 
Hie  time  increment  used  was  .05  and  NPRINT  was  5 
throughout  the  run  which  lasted  7.5  seconds. 


A  new  maneuver  was  then  set  up,  identical  to  the 
first  except  that  the  time  card,  CARD  201,  now 
contained  0.,  .05,  3.5,  .005,  3.75,  3.75  as  the 
consecutive  inputs  instead  of  0.,  .05,  7.5,  blank, 
blank,  blank  which  were  used  on  the  previous  run. 
NPRINT  on  CARD  01  was  changed  from  5  to  70.  An 
additional  CARD  211  was  input  which  had  a  J  of  31. 

Hie  number  3.5  was  in  Columns  11-20,  the  number  1.0 
in  Columns  21-30,  and  the  rest  of  the  card  blank. 

In  the  output  (see  Section  III  for  a  complete  expla¬ 
nation  of  all  outputs),  the  trim  page  was  followed 
by  the  maneuver  page  for  maneuver  time  of  0  second. 
Hie  next  time  point  for  which  output  was  given  was 
3.5  seconds  and  output  was  given  at  every  .005 
second  until  3.75  seconds.  The  result  was  no  output 
for  time  points  of  no  interest,  but  complete  coverage 
of  the  time  interval  of  interest. 
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C.  USER  TECHNIQUES 


Experience  in  using  the  program  will  increase  the  degree  of 
success  in  obtaining  trims,  maneuvers,  and  stability  calcula¬ 
tions  as  desired.  Difficulty  encountered  in  running  a  case 
may  be  due  to  a  number  of  things.  Some  of  those  related  to 
engineering  areas  are  discussed  in  Volume  I.  A  check  list  of 
the  most  likely  mechanical  and  numerical  trouble  soots,  with 
suggested  action,  follows: 

1  Input  Data 

If  a  data  deck  will  not  run,  check  for 

-  mispunched  data:  no  decimal,  extra  decimal,  wrong  field 

-  misplaced  card(s):  card  out  of  sequence,  missing  card, 
or  extra  card. 

If  a  data  deck  runs  but  results  are  obviously  wrong  (at  any 
output  level),  check  for 

an  extra  card 

special  inputs  left  in  from  a  previous  run 
wrong  number 

A  missing  card  is  very  easy  to  detect  since  it  usually  leads 
to  problems  which  generate  error  messages  by  the  operating 
system  of  the  computer.  Cards  out  of  sequence  are  usually 
very  difficult  to  detect  since  interchanging  any  pair  of  cards 
in  a  data  deck  can  generate  unique  symptoms.  An  extra  card 
usually  causes  divisions  by  zero  while  or  just  after  printing 
out  the  input  data  for  the  rotor  groups. 

2.  Difficulty  in  Trimming 

Failure  to  trim  usually  is  due  to  the  initial  approximation  of 
the  trim  vector  being  too  far  from  the  final  solution  trim 
vector.  If  a  successful  trim  has  been  computed  for  a  similar 
case,  the  new  case  can  frequently  be  trimmed  by  using  the 
variable  incrementing  feature  to  proceed  in  steps  from  the 
previous  case  to  the  new  case.  If  a  trimmed,  similar  case  is 
not  available,  a  hand  calculation  based  on  the  partial  deri¬ 
vative  matrix  output  can  aid  in  obtaining  a  closer  initial 
approximation.  Proceed  as  follows:  determine  the  three  largest 
force  or  moment  errors.  Next  find  the  largest  partial  deriva¬ 
tives  associated  with  some  group  of  three  elements  of  the  trim 
vector  in  the  three  equations  having  the  largest  errors.  Solve 
this  3x3  system  for  new  values  of  the  trim  vector  elements. 
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A  trim  failure  may  also  be  caused  by  one  variable  oscillating 
between  two  extremes  as  the  computations  go  from  one  itera¬ 
tion  to  the  next.  This  ’’bounce"  phenomenon  can  usually  be 
cured  by  using  the  average  value  in  the  initial  approximation 
and/or  reducing  the  correction  limit  and  damping  factor. 

3.  Maneuver  Diverges 

The  first  possibility  to  consider  is  bad  maneuver  inputs. 

If  the  data  are  good,  a  smaller  time  increment  should  be  tried. 
A  good  rule  of  thumb  is  that  the  time  increment  should  be  such 
that  a  blade  on  any  rotor  does  not  move  more  than  45°  around 
the  azimuth  during  one  time  increment.  It  is  also  possible 
that  the  system  may  be  inherently  unstable  and,  therefore,  the 
divergence  may  be  real. 

4.  Stability  Calculations 

Erroneous  values  for  stability  characteristics  may  result  from 
incorrect  choice  of  the  increment  used  in  numerical  approxi¬ 
mation  of  the  partial  derivatives.  If  the  increment  is  ex¬ 
tremely  small,  the  difference  computed  in  the  function  may  be 
only  due  to  round  off  error  or  may  be  zero.  If  the  increment 
is  large,  the  curvature  may  have  changed  signs  over  the  inter¬ 
val,  causing  large  errors. 

D.  SAMPLE  CASES 

In  Figures  11  through  14  are  shown  input  data  sets  for  the 
four  common  rotorcraft  configurations  for  which  the  flight 
simulation  program  is  used. 

Figure  11  shows  data  for  a  single- rotor  helicopter  flying  at 
120  knots  in  a  1.25g  coordinated  turn.  This  data  set  causes 
the  computer  program  to  find  a  trim  condition,  run  a  stability 
analysis,  and  terminate  the  problem.  The  trim  condition 
found  is  shown  in  Figure  20.  The  output  of  the  stability 
analysis  may  be  seen  in  Figures  24  through  28. 

In  Figure  12  is  the  input  data  set  for  a  compound  helicopter 
capable  of  flight  with  the  rotors  stopped.  This  particular 
data  trims  the  craft  with  the  rotors  running  at  full  speed  but 
lightly  loaded,  as  in  preparation  to  stop  and  fold  the  rotors. 
A  maneuver  to  stop  the  rotors  and  fold  the  main  rotor  was  run 
using  this  data  except  that  the  fuselage  inertias  were  made 
very  large  to  simulate  a  wind  tuhnel  test.  No  control  motions 
were  made  to  control  the  forces  and  moments.  Some  of  the 
results  of  this  case  are  shown  in  Figure  15. 


The  rotor  rpm,  thrust,  and  power  required  are  plotted  along 
with  the  torque  applied  by  the  rotor  brake  and  the  rotor 
blade  positions  after  the  rotor  rpm  is  zero.  The  brake  torque 
goes  to  zero  when  the  rotor  is  stopped  because  the  rotor  is 
assumed  to  be  locked  in  place  at  this  time. 

Figurel3  presents  data  for  a  tilting  proprotor  craft  flying 
in  helicopter  mode.  The  trim  condition  for  this  data  would 
make  a  good  starting  condition  for  a  maneuver  to  convert  from 
helicopter  mode  to  airplane  mode. 

Figure  14  lists  data  for  a  tilting  proprotor  craft  flying  in 
airplane  mode.  Starting  with  this  trim  condition  a  maneuver 
may  be  run  to  reconvert  to  helicopter  mode  or  to  stop  and  fold 
the  proprotors.  Some  results  of  the  second  case  may  be  seen 
in  Figure  16.  The  precone  angle  is  plotted  because  this  is 
the  means  by  which  the  TFTA  rotors  are  folded  in  the  computer 
program. 
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SECTION  III 


OUTPUT  GUIDE. 


All  possible  output  of  the  computer  program  has  been  divided 
into  eight  groups  for  the  purpose  of  explanation.  The  first 
three  groups  listed  will  always  be  printed.  The  appearance  of 
the  other  items  depends  on  the  options  exercised  by  the  pro¬ 
gram  user.  A  brief  explanation  of  each  output  group  is  given 
below.  The  output  is  explained  in  detail  on  the  following 
pages. 

1 .  Input  Data  for  Trim 

The  basic  set  of  input  data  is  printed  as  a  record  of  the 
parameters  used  for  each  case. 

2.  Trim  Iteration  Page 

If  input  XIT(5)  is  one,  a  trim  iteration  page  is  printed  for 
each  attempt  to  find  a  set  of  control  and  attitude  variables 
corresponding  to  a  trimmed  ship.  If  a  jrim  solution  is 
achieved,  most  of  this  information  is  superfluous.  However, 
if  a  trim  condition  cannot  be  found,  this  output  can  be  very 
useful  for  analyzing  the  problem  and  finding  a  better  set  of 
initial  values  for  control  settings,  fuselage  angles,  and 
flapping  angles. 

3.  Final  Trim  Page 

A  final  trim  page  is  always  printed  indicating  whether  or  not 
a  trim  condition  was  found  and  giving  the  general  conditions 
for  the  final  iteration  in  either  case. 

4.  Input  Data  for  Maneuver 

If  a  maneuver  is  called  for,  the  inputs  controlling  the  maneu¬ 
ver  are  printed. 

5.  Typical  Maneuver  Page 

At  desired  time  points  during  a  maneuver  a  number  of  variables 
such  as  positions,  velocities,  and  accelerations  are  printed. 

6.  Time  History  Plots 

Time  history  plots  of  variables  selected  by  the  user  may  be 
obtained  at  the  end  of  a  maneuver. 


7.  Output  of  Stability  Analysis  Routine 

During  the  stability  analysis  much  of  the  information  obtained 
in  computing  the  partial  derivatives  is  printed  out.  At  the 
completion  of  the  stability  analysis,  the  roots  of  the  charac¬ 
teristic  equations  with  controls  fixed  and  a  frequency  response 
transfer  function  are  printed. 

8 .  Least  Squares  Curve  Fit 

Following  a  maneuver  with  a  sinusoidal  control  motion  the 
curve  fit  routine  prints  coefficients,  phase  angles, and  cor¬ 
relation  factors  to  indicate  how  well  selected  variables 
follow  the  control  motion.  This  program  option  may  also  be 
used  to  find  the  mean  and  oscillatory  values  for  rotor  thrust, 
H-force,etc.  See  Figures  32  through  34. 

A.  SIGN  CONVENTIONS 

1 .  Cartesian  Coordinates 

Many  of  the  forces,  displacements,  velocities,  accelerations, 
moments,  rotational  velocities,  and  rotational  accelerations 
computed  in  this  program  are  expressed  in  cartesian  coordinates 
in  one  of  the  three  reference  systems  defined  below.  In  all 
cases  the  sense  of  the  moments  and  rotations  is  found  by  apply¬ 
ing  the  right  hand  rule  to  the  coordinate  axes. 

a.  Ground  Reference 

In  this  reference  system,  X  is  positive  due  north,  Y  is 
positive  to  the  east,  and  Z  is  positive  downward. 

b.  Fuselage  Reference 

This  coordinate  system  is  attached  to  the  rotorcraft  at  the 
center  of  gravity.  The  positive  X  axis  is  pointed  forward; 
the  positive  Y  axis  points  to  the  pilot’s  right;  and  th>. 
positive  Z  axis  is  pointed  toward  the  bottom  of  the  ro  orcraft. 

c.  Rotor  Shaft  Reference 
(1)  Z  Direction 

The  Z  axis  in  shaft  reference  lies  along  the  axis  of  the 
shaft.  Rotor  thrust  is  always  expressed  in  shaft  reference 
as  the  component  of  rotor  force  along  the  Z  axis.  The  positive 
direction  for  rotor  thrust  may  be  found  in  Table  IV.  The  posi¬ 
tive  Z  direction  is  in  the  opposite  direction  of  rotor  thrust. 
Rotor  induced  velocity  is  positive  in  the  direction  of  the 
positive  Z  axis. 
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(2)  Y  Direction 


The  positive  Y  axis  in  shaft  reference  is  in  the  direction  of 
positive  Y- force.  See  Table  IV,  The  Y-axis  is  perpendicular 
to  the  plane  of  the  fuselage  reference  X  and  Z  axes  for  all 
rotor  reference  systems  except  one.  For  the  tail  rotor  of  a 
single  rotor  helicopter,  the  Y  axis  is  perpendicular  to  the 
fuselage  reference  X-Y  plane. 

( 3)  X-Direction 

The  X  axis  is  parallel  to  the  fuselage  reference  X-Z  plane. 
The  positive  X  direction  is  opposite  the  positive  direction 
for  H-force  indicated  in  Table  IV. 


Table  IV. 

ROTOR  DESIGNATIONS  AND  FORGE  SIGN  CONVENTIONS 


Single 

Rotor 

Helicopter 

Tandem 

Rotor 

Helicop ter 

Prop -Rotor 
Aircraft 
(Helicopter 
Mode) 

Prop- Rotor 
Aircraf  t 
(Airplane 
Mode) 

Rotor  1  1 

Designation 

MAIN 

FORWARD 

RIGHT 

RIGHT 

Thrust 

Up 

Up 

Up 

Forward 

H-Force 

Aft 

Aft 

Aft 

Up 

Y-Force 

Right 

Right 

Right 

Right 

Rotor  2 

Designation 

TAIL 

AFT 

LEFT 

LEFT 

Thrus  t 

Right 

Up 

Up 

Forward 

H-Force 

Aft 

Aft 

Aft 

Up 

Y-Force 

Down 

Right 

Right 

Right 

Note:  Right  and  left  indicate  the  pilot1 s  right  and 
left  respectively. 


2 .  Euler  Angles 

The  Euler  angles  used  here  give  the  angular  orientation  of  the 
fuselage  coordinate  system  with  respect  to  ground  reference. 
They  are  a  set  of  three  ordered  angles  as  follows: 


Psi  is  a  right  handed  rotation  about  the  ground  reference 
Z  axis. 


Theta  is  a  right  handed  rotation  about  the  Y  axis  which 
has  been  rotated  through  psi  previously 


-  Phi  is  a  right  handed  rotation  about  the  X  axis  resulting 
from  the  two  prior  rotations. 


3.  Flapping  Angles  and  Moments 


The  flapping  angles  are  a]_s  and  bis  in  shaft  reference.  With¬ 
out  coning,  the  flapping  angle  may  be  found  in  terms  of  a-, 
and  b.  from  the  following  expression: 


Flapping 


-a,  cos  *1*  -  b.  sin  ^ 

Is  Is 


Where  ^  is  the  aximuth  location.  Zero  ^  is  due  aft  and  posi¬ 
tive  rotation  is  counterclockwise  when  viewed  from  above. 
Flapping  angle  in  general  is  positive  upward  so  fore  and  aft 
flapping,  ais  *  is  positive  when  a  blade  at  ^  =  0  is  down, 
that  is,  it  makes  an  angle  of  less  than  ninety  degrees  with  the 
rotor  shaft.  Lateral  flapping,  bis  ,  is  positive  when  a  blade 
at  ^  =  90  degrees  is  do.Ti. 


A  positive  fore  and  aft  'lapping  moment  is  one  which  follows 
the  right  hand  rule  about  a  line  at  =  270?  A  positive 
lateral  flapping  moment  is  a  right  hand  moment  about  a  line 
at  i|»  =  0°. 


4.  Control  Positions 


a .  Percent  or  Inches 

(1)  Collective  Stick 

Zero  percent  collective  is  full  down.  Positive  stick  motion 
in  percent  or  inches  is  upward . 

(2 )  Fore  and  Aft  Cyclic  Stick 

Zero  percent  fore  and  aft  is  full  aft.  Positive  stick  motion 
in  percent  or  inches  is  forward . 

(3)  Lateral  Cyclic  Stick 

Zero  percent  lateral  cyclic  is  full  left.  Positive  stick 
motion  in  percent  or  inches  is  to  the  right. 

(4)  Pedals 

Zero  percent  pedal  is  full  right.  Positive  pedal  motion  in 
percent  or  inches  is  to  the  left.  That  is,  positive  pedal 
tends  to  make  the  rotorcraft  yaw  nose  left. 


b .  Radians  or  Degrees 

When  control  positions  are  expressed  in  radians  or  degrees, 
these  values  apply  to  the  pitch  angles  measured  on  the  rotor 
itself.  A  neutral  stick  position  is  calculated  based  on  the 
range  of  the  stick,  the  range  of  the  control,  and  the  speci¬ 
fied  limit  on  the  control  with  the  stick  at  zero  percent. 

5 .  Miscellaneous  Quantities 

a .  Mast  Tilt  Angle 

Positive  mast  tilt  angle  for  all  rotor  configurations  is  a 
right  hand  rotation  about  the  negative  Y  axis  in  shaft 
reference . 


b .  Climb  Angle 

The  climb  angle  is  the  angle  of  the  flight  path  relative  to 
the  X-Y  plane  in  ground  reference.  It  is  positive  if  the 
rotorcraft  is  climbing. 

c .  Heading  Angle 

This  angle  is  the  direction  of  the  flight  path  on  the  compass. 
Zero  heading  is  due  north,  along  the  ground  reference  X  axis. 
A  heading  of  ninety  degrees  is  due  east. 

d .  Rotor  Rotation 

Rotor  number  one  always  rotates  in  a  right  hand  sense  about  a 
positive  thrust  vector.  Rotor  number  two  always  rotates  in 
the  opposite  direction. 


e .  Fuselage  Angles  of  Attack 

The  fuselage  angles  of  attack  are  based  on  the  fuselage 
reference  velocities,  and  gust  velocities. 

Pitch  angle  of  attack  =  tan”*-  (  el~oc^^  ^ 

°  \  X  velocity  / 

Yaw  angle  of  attack  =  tan'1  (fc  Velocity) 

Sideslip  angle  is  the  negative  of  yaw  angle  of  attack. 

f .  Gust  Velocities 

(1 )  Forward  Component 

The  forward  component  of  gust  velocity  is  positive  if  the 
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gust  is  moving  in  the  positive  X  direction  relative  to  the 
fuselage . 

(2)  Lateral  Component 


The  lateral  component  of  gust  velocity  is  positive  if  the 
gust  is  moving  in  the  positive  Y  direction  relative  to  the 
fuselage . 

(3)  Vertical  Component 

The  vertical  component  of  gust  velocity  is  positive  if  the 
gust  is  moving  in  the  positive  Z  direction  relative  to  the 
fuselage  . 

g .  Acceleration  Levels  in  Gs 


(1)  Forward 

Forward  acceleration  is  positive  in  the  positive  X  direction 
in  fuselage  reference. 

(2 )  Lateral 

Positive  lateral  g  level  is  to  the  pilots  left,  in  the  nega¬ 
tive  fuselage  Y  direction. 

(3)  Vertical 

Positive  vertical  g  level  is  upward,  in  the  negative  fuselage 
Z  direction.  For  straight  and  level  flight  the  vertical  g 
level  is  one  (1.00). 

B.  INPUT  DATA  FOR  TRIM 

1.  Disposition  of  Output  (Figure  17) 

Message  on  lead  card  for  name  and  location  of  program  user 
is  printed  on  lead  page  repeatedly. 

2 .  Problem  Identification 

The  value  of  the  primary  control  variable,  NPART,  the  problem 
identification  number,  IPSN,  and  the  three  cards  of  comments 
appear  at  the  beginning  of  each  problem.  See  Figure  18. 

3.  Input  Data  (Figure  18) 

The  fourteen  basic  groups  of  data  are  printed  exactly  as  they 
are  input . 


4 .  Check  Aerodynamic  Inputs 

Several  of  the  input  aerodynamic  inputs  are  changed  if  they 
do  not  have  a  reasonable  value.  An  error  message  is  printed 
explaining  the  action  taken.  See  Figure  18  following  input 
data. 

5 .  Non-Standard  Trim  Conditions 

If  the  rotorcraft  is  to  be  trimmed  in  a  coordinated  turn,  a 
pull-up,  or  a  push-over,  information  is  printed  concerning 
these  conditions  (Figure  18). 

C.  TRIM  ITERATION  PAGE  (Figure  19) 

1.  VAR  (i) 

This  row  gives  the  current  value  of  the  ten  variables  which 
are  changed  in  order  to  trim  the  rotorcraft. 

VAR(l)  =  Collective  Stick  Position  in  Percent 

VAR(2 )  =  Fore  and  Aft  Cyclic  Stick  Position  in  Percent 

VAR(3 )  =  Lateral  Cyclic  Stick  Position  in  Percent 

VAR(4)  =  Pedal  Position  in  Percent 


VAR(5)  =  Fuselage  Pitch  Angle  in  Degrees 

VAR(6)  =  Fuselage  Roll  or  Yaw  Angle  in  Degrees  [choice 
of  roll  or  yaw  is  made  by  user,  XIT(15)] 

VAR(7 )  =  Main,  forward,  or  Right  Rotor  Fore  and  Aft 
Flapping  Angle  in  Degrees  (a^s) 

VAR(8)  =  Main  Forward,  or  Right  Rotor  Lateral  Flapping 
Angle  in  Degrees  (b]^s ) 

VAR(9 )  =  Tail  Aft,  or  Left  Rotor  Fore  and  Aft  Flapping 
Angle  in  Degrees  (a^s ) 

VAR(10)=  Tailt  Aft,  or  Left  Rotor  Lateral  Flapping  Angle 
in  Degrees  (b^fi) 

2 .  Rotor  Performance 

These  two  rows  give  the  following  quantities  for  the  two 
rotors : 


-  Thrust  in  shaft  reference  (lb); 

-  H-force  in  shaft  reference  (lb); 

-  Y-force  in  shaft  reference  (lb); 

-  Torque  in  shaft  reference  (ft-lb) ; 

-  Average  Induced  Velocity  (ft/sec) ; 

-  Coning  Angle  (deg); 

-  Also  thrust  provided  by  the  jets  if  present. 

3.  Force  and  Moment  Summary 

This  block  of  output  shows  the  contribution  of  each  element  of 
the  rotorcraft  to  the  total  forces  and  moments.  The  X -force, 
Y-force,  Z-force,  Roll  Moment,  Pitch  Moment,  and  Yaw  Moment 
are  in  fuselage  reference  with  the  forces  in  pounds  and  the 
moments  in  foot-pounds.  The  rotor  flapping  moments  are  in 
shaft  reference  with  the  units  being  foot-pounds.  The  columns 
of  forces  and  moments  are  presented  as  labeled  in  the  following 
order: 

-  Total  caused  by  all  elements 

-  Right  Wing 

-  Left  Wing 

-  Elevator  (Horizontal  Stabilizer) 

-  Fuselage 

-  Right  or  Center  Jet 

-  Left  Jet 

-  Main,  Forward,  or  Right  Rotor 

-  Tail,  Aft,  or  Left  Rotor 

-  Gun  (Non-zero  for  a  Maneuver  only) 

-  Fin  (Fin/Rudder) 

-  Weight/^ain  or  \\yAt  P  .tor  Torque  (QMR) 

The  forces  m  this  column  are  caused  by  the  weight; 
the  moments  <ire  caused  by  main,  forward,  or  right 
rotor  torque . 
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-  Tail,  aft,  or  left  rotor  torque  (QTR) . 

4.  Partial  Derivative  Matrix 

This  matrix  gives  the  partial  derivative  of  each  force  and 
moment  with  respect  to  each  of  the  iteration  variables  . 

The  units  on  the  force  derivatives  are  pounds  per  radian, 
and  on  the  moment  derivatives  foot-pounds  per  radian.  For 
the  controls  the  angles  are  measured  on  the  rotor.  If 
input  XIT(3)  =  0  this  matrix  is  computed  and  printed  at 
every  iteration,  otherwise  it  appears  every  fifth  iteration. 

5 .  Corrections 

The  line  labeled  "CORRECTIONS'*  gives  the  computed  changes 
in  the  iteration  variables  in  radians.  They  are  in  the 
same  order  as  the  VAR(I)  and  the  partial  derivative  rows 
If  any  of  these  computed  correct: ions  are  greater  than  the 
maximum  allowed  by  the  inputs  in  the  iteration  limits  group, 
the  entire  row  is  multiplied  by  a  ratio  which  will  make  all 
corrections  within  the  allowable  range  and  this  ratio  is 
printed  along  with  the  number  of  the  iteration  variable  which 
determined  it.  The  corrections  are  then  added  to  the  itera¬ 
tion  variables  to  determine-  the  values  for  the  next  iteration. 
It  should  be  noted  again  that  the  "CORRECTIONS"  are  in  radians 
and  not  the  same  units  as  the  VAR(I)  printed. 

D.  FINAL  TRIM  PAGE  (Figure  20) 

1 .  Problem  Identification 

This  is  the  same  as  mentioned  in  B.  "INPUT  DATA  FOR  TRIM.” 

2 .  Trim  Condition 

A  one-line  message  is  printed  indicating  whether  or  not  the 
rot or craft  is  in  a  stable  condition. 

j .  Rotor  Conditions 

At  the  upper  left  of  this  page  are  the  basic  rotor  performance 
quantities.  Most  of  this  information  is  clearly  marked  and 
self-explanatory.  The  collective  pitch  is  determined  at  the 
root  of  the  blade.  The  flapping  angles,  thrust,  H-force, 
and  Y-force  are  in  shaft  reference. 

4 .  Jet  Thrust 

Below  the  rotor  data  are  the  values  of  jet  thrust.  Jet  thrust 
may  be  varied  with  the  collective  stick  during  the  trim  proce¬ 
dure,  so  it  is  necessary  to  have  the  output  on  this  page. 
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5.  Percent  Control  Used 


Below  the  jet  data  the  control  positions  are  given  in  percent. 
For  the  collective  stick,  zero  percent  is  full  down.  For  the 
fore-and-aft  cyclic  stick,  zero  percent  is  full  aft.  For  the 
lateral  cyclic  stick,  zero  percent  is  full  left.  For  the 
pedals,  zero  percent  is  full  right. 

6.  Aerodynamic  Surfaces 

The  angles  of  attack  and  force  components  in  fuselage  reference 
for  the  wing,  elevator  and  fin  are  printed  in  the  upper  right 
portion  of  the  trim  page. 

7.  Fuselage  Data 

Below  the  data  on  the  aerodynamic  surfaces  is  given  the  basic 
information  on  the  aircraft  flight  condition  and  attitude. 

The  order  of  the  Euler  angle  rotations  from  fixed  reference  to 
fuselage  reference  is  yaw,  pitch,  and  roll. 

8.  Iteration  Count  and  Time 

At  the  bottom  of  the  page  the  number  of  iterations  required 
to  trim  the  rotorcraf t*  and  the  elapsed  computing  time  are 
printed. 

E„  INPUT  DATA  FOR  MANEUVER  (Figure  21) 

The  program  prints  the  contents  of  the  time  card  and  all 
maneuver  control  cards  as  they  are  input  before  the  start  of 
the  actual  maneuver.  This  serves  as  a  record  of  what  type  of 
maneuver  has  been  run  as  well  as  a  quick  way  to  check  the 
input  data. 

F.  TYPICAL  MANEUVER  PAGE  (Figure  22.) 

1.  Time 

At  the  top  of  the  page  for  each  maneuver  time  point  is  printed 
the  current  time  in  the  maneuver  and  the  total  time  used  by 
the  computer  to  this  point. 

2.  Rotor  Variables 

Most  of  the  rotor  data  of  interest  are  printed  for  each  rotor 
separately.  The  main,  forward,  or  right  rotor  data  are  printed 
first  followed  by  the  tail,  aft,  or  left  rotor  data.  These 
data  are  broken  up  into  six  major  groups. 


a.  Tip  Path  Plane  Location  and  Velocity 

The  fore  and  aft  flapping  velocity,  Q,  and  position,  A1 ;  and 
the  lateral  flapping  velocity,  P,  and  position,  Bl,  are 
printed  at  the  upper  left.  The  units  are  degrees  per  second 
and  degrees.  All  four  quantities  are  in  shaft  reference. 

b.  Hub  Velocities 

To  the  right  of  group  (1)  are  U,  V,  and  W,  which  are  the  X,  Y, 
and  Z  velocity  components,  respectively,  of  the  hub  in  shaft 
reference.  The  units  are  feet  per  second. 

c.  Conditions  for  Rotor  Blade  Number  One 

To  the  right  of  group  (2)  is  printed  information  on  the  loca¬ 
tion,  velocity,  and  acceleration  of  blade  number  one  for  the 
rotor  concerned.  The  column  headed  "PSI"  gives  the  rotational 
acceleration,  velocity,  and  location;  the  column  headed  "BEtfA" 
gives  the  flapping  acceleration,  velocity,  and  location  rela¬ 
tive  to  the  snaftT  The  units  are  consistent  in  degrees  and 
seconds . 

d .  Rotor  Forces 

To  the  right  of  group  (3)  the  forces  generated  by  the  rotor 
are  printed.  The  instantaneous  values  of  thrust,  H-force, 
and  Y-force  in  pounds  are  in  shaft  reference. 

e .  Control  Positions 

Below  group  (1)  are  given  the  collective,  fore  and  aft  cyclic, 
and  lateral  cyclic  pitch  settings  in  degrees  on  the  rotor. 

The  row  labeled  "CONTROLS"  is  determined  by  the  stick  posi¬ 
tions.  The  row  labeled  "OTHER"  is  control  contributions  caused 
by  the  equations  representing  optional  devices.  For  collective 
pitch  the  "OTHER"  control  control  contributions  may  be  caused 
by  pitch-cone  coupling,  a  prop-rotor  collective  governor,  or 
the  pitch  A.S.E.  (Automatic  Stabilization  Equipment).  For 
fore  and  aft  or  lateral  cyclic  "OTHER"  control  contributions 
may  be  caused  by  the  focused  pylon  effect  or  the  operation  of 
a  flat  tracker.  The  "TOTAL"  row  gives  the  algebraic  sum  of 
the  preceding  two  rows  for  convenience. 

f .  Miscellaneous  Data 

To  the  right  of  group  (5)  are  several  variables  of  interest 
which  cannot  be  classified  as  a  group.  The  rotor  coning 
angle  in  degrees  includes  built-in  precone  as  well  as  elastic 
effects.  "IND.V"  indicates  the  average  induced  velocity  from 
modified  momentum  theory  in  feet  per  second.  ,rM.TILT" 
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indicates  mast  tilt  angle  in  degrees.  For  a  main  rotor,  or 
prop-rotor,  or  the  aft  rotor  on  a  tandem,  configuration,  zero 
mast  tilt  is  vertical.  For  the  tail  rotor  on  the  usual  heli¬ 
copter  configuration,  zero  mast  tilt  is  to  the  pilot's  right, 
that  is,  in  the  positive  Y  direction.  "RPM"  indicates  the 
instantaneous  rotor  speed  in  revolutions  per  minute.  This 
value  will  remain  constant  unless  the  power  available  is  ex¬ 
ceeded.  The  value  of  torque  given  beside  "TORQ"  is  the  airload 
torque  on  this  rotor  in  foot-pounds.  The  horsepower  required 
to  drive  this  rotor  at  this  instant  is  indicated  by  “HP." 

The  present  position  of  the  flapping  stops  is  given  under 
"FLAP.LIM"  in  degrees.  The  stops  are  set  so  that  when  one 
blade  hits  the  upper  limit,  a  blade  180  degrees  around  the 
azimuth  will  hit  the  lower  limit  at  the  same  time.  The  cur¬ 
rent  values  of  the  flapping  springs,  which  may  vary  with 
rotor  speed,  are  given  under  "HUB  SPRINGS"  in  foot-pounds  per 
degree . 

3.  Ground  Reference  Variables 

The  position  and  velocity  of  the  rotorcraft  relative  to  the 
ground  is  given  in  several  ways. 

a.  Cartesian  Coordinates 

The  X,  Y,  and  Z  components  of  velocity  and  location  are 
printed  on  the  left  of  the  page.  The  units  are  feet  per 
second  and  feet  respectively. 

b.  Position 

The  position  of  the  craft  is  partially  defined  again  by  the 
distance  flown  and  present  altitude.  The  units  of  these 
quantities  are  feet. 

c .  Speed 

The  total  airspeed  and  net  ground  speed  are  also  given  in 
knots . 

d .  Flight  Path  Angles 

The  flight  path  heading  relative  to  due  north  is  given  in 
degrees.  The  climb  angle  relative  to  the  ground  plane  is 
given  in  degrees. 

4.  Fuselage  Reference  Variables 

The  accelerations  and  velocities  of  the  rotorcraft  are  given 
in  this  group  along  with  the  Euler  angle  locations  and  velo¬ 
cities. 
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a.  Translational  Velocities  and  Accelerations 


In  fuselage  reference  U,  V,  and  W  are  the  X,  Y,  and  Z  compo¬ 
nents  of  velocity.  The  fuselage  referenced  velocities  and 
accelerations  are  printed  under  these  headings  to  distinguish 
them  from  ground  reference  quantities.  The  units  are  consis¬ 
tent  in  feet  and  seconds. 

b.  Rotational  Velocities  and  Accelerations 

In  fuselage  reference  P,  Q,  and  R  are  the  roll  rate,  pitch 
rate,  and  yaw  rate,  respectively.  Under  the  P,  Q,  and  R 
headings  are  printed  the  angular  accelerations  and  the  angular 
velocities  in  degree  and  second  units. 

c.  Bobweight  Acceleration  and  Velocity 

The  differential  equation  for  the  bobweight  is  written  in  terms 
of  degrees  of  collective  pitch  on  the  main  rotor.  The  accel¬ 
eration  and  velocity  printed  under  "BOBWT11  are  in  these  same 
units. 

d.  Euler  Angles  from  Ground 

The  Euler  angles  given  are  'those  to  get  from  ground  reference 
to  fuselage  reference.  The  order  of  the  rotations  is  first 
yaw,  "PSI,"  (positive  nose  right),  then  pitch,  "THETA" 

(positive  nose  up),  and  then  roll,  "PHI"  (positive  left  side 
up,  right  side  down).  The  Euler  angle  velocities  are  not  the 
same  as  the  corresponding  angular  rates  R,  Q,  and  P  in  general. 
The  units  are  degrees  per  second  and  degrees. 

5.  Aerodynamic  Surfaces  and  Miscellaneous 

Below  the  fuselage  reference  variables  are  several  groups  of 
unrelated  data. 

a .  Control  Stick  Positions 

The  stick  positions  in  percent  are  printed  to  show  the  opera¬ 
tion  of  any  maneuver  control  input? 

b.  Aerodynamic  Surfaces 

The  angle  of  attack  in  degrees,  "ATK,"  and  the  lift,  "CL," 
and  drag,  "CD,"  coefficients  are  printed  for  each  of  the  aero¬ 
dynamic  surfaces,  (left  wing,  right  wing,  elevator,  and  fin/ 
rudder).  A  positive  angle  of  attack  for  a  wing  or  elevator 
will  produce  a  lift  force  which  is  upward.  A  positive  angle 
of  attack  for  the  fin/rudder  will  produce  a  force  to  the 


c .  Fuselage  Angles  of  Attack 

The  "ATKY"  label  indicates  the  yaw  angle  of  attack  of  the 
fuselage  in  degrees.  This  is  the  negative  of  the  sideslip 
angle.  The  "ATKP"  label  indicates  the  pitch  angle  of  attack 
of  the  fuselage  in  degrees. 

d.  Center  of  Gravity  Location 

The  station  iine,  buttline,  and  waterline  location  of  the 
rotorcraft  center  of  gravity  is  given  in  inches.  These  values 
may  change  as  the  rotors  are  tilted  or  folded. 

e .  Gust  Velocities 

The  forward,  lateral,  and  vertical  components  of  gust  velocity 
in  fuselage  reference  at  the  center  of  gravity  are  printed  in 
feet  per  second. 

f .  Acceleration  Levels 

The  acceleration  level  felt  by  a  sensor  at  the  center  of  gra¬ 
vity  is  printed  in  gs.  The  three  components  of  acceleration 
are  given  in  fuselage  reference . 

6.  Power  and  Related  Data 

a.  Jet  Thrust 

The  thrust  due  to  each  jet  engine  is  given  in  pounds. 

b.  Engine  Data 

Engine  torque  supplied  and  shaft  horsepower  output  are  printed 
together.  ‘This  power  is  either  the  power  required  or  the 
maximum  power  available  whichever  is  lower.  The  torque  is 
found  directly  from  this  power  value.  The  units  on  torque  are 
foot-pounds . 

c.  Total  Horsepower  and  Rotor  Brake 

The  rotor  brake  torque  is  given  in  foot-pounds.  This  value 
is  computed  so  that  one  blade  of  the  rotor  will  stop  at  the 
position  specified  on  the  maneuver  input.  A  negative  value 
u  iy  be  computed  in  the  last  instant  before  the  rotor  stops 
but  not  for  more  than  one  time  point.  The  total  horsepower 
required  is  obtained  from  the  total  torque  due  to  the  rotors 
and  the  rotor  brake.  A  difference  between  this  quantity  and 
the  engine  shaft  horsepower  indicates  that  available  power 
has  been  exceeded. 
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7.  Force  and  Moment  Summary 

This  block  of  output  is  as  described  in  C.  TRIM  ITERATION  PAGE. 

8.  Instantaneous  Rotor  Flapping 

For  each  blade  on  each  rotor  the  azimuth  location,  the  local 
flapping  acceleration,  the  local  flapping  velocity,  and  the 
local  flapping  angle  are  printed  in  tabular  form.  See  Figure 
22.  The  units  are  consistent  in  degrees  and  seconds. 

G.  TIME  HISTORY  PLOTS  (Figure  23) 

1 .  Problem  Identification 

This  printout  is  as  described  in  B.  INPUT  DATA  FOR  TRIM. 

2 •  Variable  Ranges  and  Scales 

The  lower  and  upper  limits  on  the  plot  scale  are  given  for 
each  of  the  three  variables  which  may  be  plotted  by  each  plot 
card.  The  scale  in  units  per  inch  is  also  g'.ven. 

3.  Variables  Plotted 

The  plot  symbols  used  are  the  numbers  1,  2,  and  4.  The  vari¬ 
able  corresponding  to  each  symbol  and  its  units  are  printed 
as  part  of  the  plot  heading. 

4.  General  Comments 

The  plot  symbols  1,  2,  and  4  are  used  so  that  every  plot  point 
can  be  identified.  If  more  than  one  symbol  falls  on  the  same 
print  position,  the  symbol  printed  is  the  sum  of  these.  An 
explanation  of  the  print  symbols  3,  5,  6,  and  7  is  included 
as  part  of  each  plot  heading. 

A  good  deal  of  judgment  is  required  to  select  the  best  plot 
ranges.  If  the  selected  range  is  too  large,  small  variations 
of  the  plotted  variable  will  be  hidden.  If  the  selected  range 
is  too  small,  many  points  may  be  lost  off-scale,  and  overall 
trends  may  be  hidden.  The  range  may  also  be  shifted  to  the 
high  or  low  side  so  that  no  points  at  all  are  plotted  in  ex¬ 
treme  cases.  Any  points  which  do  not  fall  within  the  ten- 
inch  range  allowed  on  the  plot  page  are  not  printed. 

The  maneuver  time  is  printed  along  the  side  of  the  plot  so 
that  the  full  maneuver  output  may  be  easily  located  for  any 
given  plot  point.  If  the  time  increment  is  changed  during 
a  maneuver,  there  will  be  an  apparent  change  in  the  time  scale 
at  this  point. 
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Each  plot  card  is  independent  of  all  other  plot  cards.  Thus, 
if  desired,  one  variable  may  be  plotted  on  more  than  one  plot. 
One  example  which  has  proved  useful  is  rotor  azimuth  position. 

The  solitary  number  following  the  plot  is  the  elapsed  time 
in  the  computer  to  make  this  plot. 

H.  OUTPUT  OF  STABILITY  ANALYSIS  ROUTINE 

The  operation  of  the  stability  analysis  depends  on  the  numer¬ 
ical  evaluation  of  a  number  of  partial  derivatives.  The 
partial  derivatives  appear  in  the  equations  of  motion  for  the 
rotorcraft.  A  frequency  analysis  is  made  on  the  equations  of 
motion  with  controls  fixed  and  following  impulse  inputs  to  the 
controls.  "S’*  is  the  Laplace  operator  as  used  here  as  well  as 
in  the  computer  output. 

I.  Control  and  Attitude  Partial  Derivatives  (Figure  24) 

A  7  by  6  matrix  similar  to  the  partial  derivative  matrix 
printed  during  the  trim  procedure  is  printed.  The  first  four 
rows  give  the  partial  derivatives  of  the  six  forces  and  mo¬ 
ments  (X -FORCE,  Y-FORCE,  Z-FORCE,  YAW  MOM,  PITCH  MOM,  and  ROLL 
MOM),  in  fuselage  reference,  with  respect  to  each  of  the  four 
control  motions.  These  quantities  are  in  pounds  per  inch  of 
control  and  foot-pounds  per  inch  of  control.  The  fifth  and 
sixth  rows  printed  are  the  partial  derivatives  of  the  six 
forces  and  moments  with  respect  to  pitch  angle  and  roll  angle 
respectively.  These  quantities  are  in  pounds  per  radian  and 
foot-pounds  per  radian.  The  bottom  row,  labeled  ”  -  ERROR," 
is  printed  because  of  a  programming  simplification  and  should 
be  disregarded. 

2.  Velocity  Partial  Derivatives  (Figure  25) 

The  next  six  pages  of  output  contain  detailed  information  used 
for  the  calculation  of  the  partial  derivatives  with  respect  to 
the  six  fuselage  velocity  components.  The  partial  derivatives 
are  taken  in  the  same  order  as  the  variables  (VAR(I)  which  are 
listed  below. 

a.  VAR(I) 


In  the  stability  analysis  output  the  VAR(I)  are  the  six  fuse¬ 
lage  reference  velocity  components  in  the  following  order: 

VAR(l)  =  U,  velocity  in  the  X  direction  (feet  per 
second) 

VAR(2)  =  W,  velocity  in  the  Z  direction  (feet  per 
second) 


. 


VAR(3)  =  Q,  pitch  rate  (degrees  per  second) 

VAR(4)  =  V,  velocity  in  the  Y  direction  (feet  per 
second) 

VAR(5)  =  P,  roll  rate  (degrees  per  second) 

VAR(6)  =  R,  yaw  rate  (degrees  per  second) 

b.  Rotor  Performance 

These  two  rows  are  as  described  in  part  2  under  section  C 
TRIM  ITERATION  PAGE. 

c.  Force  and  Moment  Summar 

This  block  of  output  is  as  described  in  part  3  under  section 
C.  TRIM  ITERATION  PAGE.  The  forces  and  moments  printed  here 
are  computed  after  the  small  increment  in  the  pertinent  ve¬ 
locity  component  has  been  made.  The  rotor  flapping  angles 
have  also  been  adjusted  so  that  the  flapping  moments  are  the 
same  as  before  the  velocity  increment. 

d.  Delta/Force  and  Moment  Summar 


This  block  of  output  presents  the  changes  in  the  force  and 
moment  contributions  in  exactly  the  same  format  as  the  full 
FORCE  AND  MOMENT  SUMMARY.  Any  number  in  this  block  is  ob¬ 
tained  by  taking  the  corresponding  value  from  the  FORCE  AND 
MOMENT  SUMMARY  immediately  above  less  the  corresponding  value 
at  the  trim  condition  or  at  the  current  maneuver  time  point. 

If  the  stability  analysis  is  called  during  a  maneuver,  the 
forces  and  moments  are  recomputed  based  on  the  quasi-static 
rotor  analysis.  Thus  the  base  values  for  the  partial  deriva¬ 
tives  of  the  rotor  forces  will  probably  not  be  the  same  as 
those  printed  on  the  maneuver  page.  The  maneuver  printout 
contains  instantaneous  values  whereas  the  stability  analysis 
is  based  on  mean  values.  Note  that  the  rotor  flapping  moments 
remain  unchanged  from  the  previous  FORCE  AND  MOMENT  SUMMARY 
printout.  This  is  because  the  rotor  flapping  moments  are 
allowed  to  return  to  the  input  values  by  adjusting  the  flap¬ 
ping  angles  following  the  small  velocity  change. 

e.  Normalized  Delta/Force  and  Moment  Summa 


This  block  of  data  presents  basically  the  same  information  as 
the  block  above.  However,  the  values  are  all  normalized  so 
that  the  total  changes  in  the  forces  and  moments  are  one  hun¬ 
dred.  Thus  the  quantities  printed  represent  a  percent  of  the 
total  changes.  Again  the  rotor  flapping  moments  are  not 
changed • 


3.  Stability  Partial  Derivative  Matrices  (Figure  26) 


On  this  page  is  printed  a  summary  of  the  partial  derivatives 
computed  from  the  data  on  the  six  previous  pages.  Each  row 
gives  the  partial  derivatives  of  some  force,  moment,  or  flap¬ 
ping  angle  as  labeled,  with  respect  to  the  linear  and  angular 
velocities,  U,  W,  Q,  V,  P,  and  R.  The  units  are  consistent 
in  feet,  pounds,  radians,  and  seconds. 

4.  Longitudinal  Mode  (Figure  27) 

The  longitudinal  mode  deals  with  forward  velocity,  U,  fuselage 
angle  of  attack,  ALPHA,  and  fuselage  pitch  angle,  THETA,  as 
dependent  variables.  The  variables  which  act  as  forcing  func¬ 
tions  are  the  fore-and-aft  cyclic  stick  and  the  collective 
stick. 

a.  Equations  of  Motion 

The  coefficients  for  the  linearized  differential  equations  of 
motion  for  the  summation  of  X-forces,.  the  summation  of  Z- 
forces,  and  the  summation  of  pitching  moments  are  printed  in 
both  dimensional  and  nondimensional  form.  The  column  headings 
indicate  the  dependent^  variable  and  the  order  of  the  derivative. 
In  the  Laplace  operator  notation  used,  S  indicates  the  first 
derivative  and  S**2  indicates  the  second  derivative.  The 
columns  headed  "F/A  Cyclic"  and  "COLLEC."  indicate  the  effec¬ 
tiveness  of  the  controls,  and  the  signs  used  assume  they  are 
on  the  right-hand  side  of  the  equations.  That  is,  the  "=" 
sign  should  be  placed  between  the  "THETA"  column  and  the  "F/A 
CYCLIC"  column. 

b.  Controls  Fixed  Analysis 

The  characteristic  equations  are  expanded  into  polynominal 
form  and  factored  to  give  complex  roots,  S  =  a  +  ib.  The 
values  of  a  and  b  are  printed  under  "REAL"  and  "IMAG."  re¬ 
spectively.  Only  the  positive  imaginary  part  of  complex 
roots  is  printed. 

The  ,fPERIOD"  and  "FREQUENCY"  are  non-zero  only  for  roots  with 
a  non-zero  imaginary  part  because  they  have  no  meaning  other¬ 
wise.  These  quantities  are  given  in  seconds  and  cycles  per 
second  respectively. 

The  column  headed  "T  *  HALF  -  DBL"  depends  only  on  the  value 
of  the  real  part  of  the  root.  If  the  real  part  is  negative, 
the  time  to  half  amplitude,  in  seconds,  is  printed.  If  the 
real  part  is  positive,  the  time  to  double  amplitude,  in  sec¬ 
onds,  is  printed. 
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The  mode  shape  for  each  of  the  roots  is  presented  under 

"FWD . VEL/THETA"  and  "ALPHA/THETA. "  The  output  lists  a  relative 

magnitude,  "MAGN"  and  a  phase  angle  difference,  "PHASE." 

The  column  labeled  "LEAD  COEF."  presents  the  coefficient  of 
the  highest  power  of  S  in  the  polynominal  expansion  of  the 
characteristic  equations. 

c .  Frequency  Response  Roots 

The  frequency  response  roots  for  the  six  combinations  of  de¬ 
pendent  variable  and  independent  variable  are  found.  The 
roots  may  be  listed  in  the  following  form: 

S  =  REAL1  +  i  I MAGI 

S  =  REAL 2  +  i  IMAG2 

S  =  REAL3  +  i  IMAG3 

If  one  of  the  imaginary  parts  is  non-zero,  the  complex  con¬ 
jugate  root  will  be  printed. 

The  value  printed  under  "GAIN"  is  the  constant  term  in  the 
frequency  response  polynominal . 

d.  Numerators  for  Transfer  Functions 

The  numerator  for  the  Laplace  transfer  function  is  printed 
for  each  of  the  six  combinations  of  dependent  variable  and 
independent  variable.  The  static  gain  and  three  quadratic 
factors  may  be  printed.  If  both  the  "TAU"  and  the  "DAMP" 
coefficients  are  zero  for  one  factor,  then  that  factor  may 
be  omitted  from  the  numerator.  The  numerator  of  the  Laplace 
transfer  function,  N(S),  may  be  obtained  from  the  output  as 
follows : 

N(S)  =  (STATIC  GAIN)  *  (TAU  S2  +  DAMP  S  +  1) 

*( TAU  S2  +  DAMP  S  +  1)  *  (TAU  S2  +  DAMP  S  +  1) 


e.  Denominator  to  be  used  for  the  six  transfer  functions  is 
derived  from  the  roots  of  the  characteristic  equations  with 
the  controls  fixed.  It  is  printed  in  the  same  format  as  the 
numerators.  The  value  under  "STATIC  GAIN"  in  this  case 
should  be  disregarded.  The  denominator  of  the  Laplace  trans¬ 
fer  function,  D(S),  may  be  obtained  from  the  output  as  fol¬ 
lows  : 

D(S)  =  (TAU  S2  +  DAMP  S  +  1)  *  (TAU  S2  +  D/JiP  S  +  1) 

*  (TAU  S2  +  DAMP  S  +  1) 


A  complete  transfer  function  is  formed  by 

f .  General  Information 

At  the  bottom  of  the  page  are  several  notes  concerning  the 
units  used  in  the  output  for  quick  reference.  The  "NATURAL 
UNITS"  referred  to  are  either  feet  per  second  or  radians. 

The  last  line  gives  a  value  of  T* ,  the  nondimensionalizing 
factor  for  time.  For  the  longitudinal  mode,  T*  equals  the 
mean  aerodynamic  chord  divided  by  twice  the  total  velocity. 

5.  Lateral  Mode  (Figure  28) 

The  lateral  mode  deals  with  side.' lip  angle,  BETA,  roll  angle, 
PHI,  and  yaw  rate,  R,  as  dependent  variables.  The  variables 
which  act  as  forcing  functions  art  lateral  cyclic  stick  and 
the  pedals. 

Note  that  yaw  angle  and  yaw  rate  are  positive  nose  right 
while  sideslip  angle  is  positive  nose  left. 

The  output  for  the  lateral  mode  is  in  a  form  identical  to 
that  of  the  longitudinal  mode.  Only  the  differences  will  be 
mentioned  below. 

a.  Equations  of  Motion 

The  equations  of  motion  for  the  lateral  mode  are  for  the 
summation  of  Y-forces,  the  summation  of  rolling  moments,  and 
the  summation  of  yawing  moments  in  that  order. 

b.  Static  Gain 

The  units  on  the  static  gains  are  radians  per  inch  of  stick 
or  radians-per-second  per  inch  of  stick. 

c.  T* 

For  the  lateral  mode,  T*  is  the  wing  span  divided  by  twice  the 
total  velocity. 

I.  VECTOR  ANALYSIS  DATA 

This  program  option  gives  a  harmonic  analysis  of  selected 
variables  following  the  completion  of  a  maneuver. 


1.  Problem  Identification 

This  output  is  as  described  in  B.  INPUT  DATA  FOR  TRIM. 


The  maneuver  time  at  which  the  curve  fit  starts  is  given.  All 
time  points  prior  to  this  time  are  disregarded  by  the  curve 
fit  procedure.  The  frequency  used  in  the  curve  fit,  OMEGA,  is 
given  in  cycles  per  second.  The  curve  fit  function,  F(T),  is 
expressed  in  general  form: 

F(T)  =  AMPLITUDE  *  SIN( OMEGA  *  T  +  PHASE  ANGLE)  + 

CONSTANT 

(where  T  is  time  as  measured  during  the  maneuver). 

3.  Variable ,  Amplitude ,  Phase  Angle,  and  Constant  ( Figures 
29  and  32) 

Below  the  general  equation  are  five  columns  as  follows: 

a.  Variable 

In  this  column  the  variable  being  curve  fit  is  identified, 
and  its  units  are  given. 

b.  Amplitude 

This  number  may  be  substituted  into  the  general  equation  for 
AMPLITUDE.  The  units  are  those  given  under  VARIABLE. 

c.  Phase  Angle 

This  number  may  be  substituted  into  the  general  equation  for 
PHASE  ANGLE.  The  units  are  degrees  as  labeled. 

d.  Constant 

This  quantity  may  also  be  substituted  directly  into  the  gene¬ 
ral  equation.  The  units  are  those  given  under  VARIABLE. 

e .  Coef  of  Corr 

This  denotes  coefficient  of  correlation  and  is  a  measure  of 
how  well  the  variable  considered  is  fit  by  a  sinusoidal  varia¬ 
tion  at  the  frequency  selected.  A  number  greater  than  0.95  in 
this  column  indicates  a  reasonably  good  fit.  A  number  smaller 
is  generally  caused  by  other  frequency  content  or  transient 
conditions. 

4.  Problem  Identification 

At  the  top  of  the  following  page  the  problem  identification  is 
repeated. 
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5.  Amplitude  and  Phase  Angle  Comparisons  (Figures  30  and  33) 

A  comparison  of  the  magnitudes  and  phase  angles  between  vari¬ 
able  vectors  may  be  made  for  selected  pairs  of  variables.  The 
variables  compared  are  labeled  as  variable  A/variable  B.  The 
variable  identifications  used  are  the  same  as  those  used  on  the 
previous  page  and  for  the  plot  headings.  The  amplitude  ratio 
printed  is  AMPLITUDE  A  divided  by  AMPLITUDE  B.  The  phase  angle 
difference  is  PHASE  ANGLE  A  minus  PHASE  ANGLE  B. 

6.  Problem  Identification 

Following  the  amplitude  and  phase  angle  comparisons,  the  pro¬ 
gram  skips  to  the  top  of  the  next  page  and  again  prints  the 
problem  identification  heading. 

7.  Variable  "A”  as  a  Linear  Combination  of  Variables  "B"  and 
"C"  (Figures  3l  and  3*0 

If  all  the  selected  variables  are  viewed  as  vectors  rotating 
at  the  same  rotational  speed,  OMEGA,  any  one  variable  may  be 
expressed  as  a  linear  combination  of  two  other  variables  and 
a  constant  as  long  as  the  phase  angle  between  the  two  variables 
is  not  0  or  180  degrees.  This  relationship  is  given  generally 
in  the  heading  as  "A  =*KB  *  B  +  K0  *  C  +  KD. " 

Here  A,  B,  and  C  are  the  variables  concerned.  The  variable 
identification  phrase  is  printed  for  each  in  the  output.  KB, 
KC,  and  KD  are  constants  determined  by  the  program  and  printed 
in  the  column  labeled  COEFFICIENT.  In  this  row  for  variable 
B,  the  coefficient  is  KB;  in  the  row  for  variable  C  the  coef¬ 
ficient  is  KC;  and  in  the  unlabeled  row,  which  has  CONSTANT, 
the  coefficient  is  KD. 

8.  Time  Used 

At  the  completion  of  the  vector  analysis  routine  the  time  used 
in  the  vector  analysis  process  is  printed  along  with  the  total 
elapsed  computing  time. 

J.  ROTOR  AIRLOAD  DATA 

At  maneuver  time  points  specified  in  the  airload  times  group, 
the  airload  distribution  on  the  rotors  is  printed.  The  option 
to  print  airloads  for  either  or  both  rotors  is  determined  by 
the  input  value  of  NVARC  on  data  card  01.  A  set  of  airload 
data  is  printed  for  each  blade  on  the  rotor  or  rotors  selected. 
See  Figure  35. 

1.  Variables  Which  Are  Independent  of  Radius 
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a.  "AZIMUTH" 

degrees?Uth  P°Siti°n  °f  the  blade  in  Option  i8  given  in 

b.  "U-SHAFT" 

givenh?ntfntepeerCLcond”POnent  °f  the  r°t0r  hub  vel°^ty  is 

c.  1 V- SHAFT" 

givenh?ntf«fepeerC:econd?POI’ent  °f  ^  r°t0r  hub  vel-ity  is 

d-  "W- SHAFT" 

givenh!ntfntep:rCsecond?P°nent  °f  the  r0t°r  hub  velocity  is 
e.  "XK" 

de termines^the ^variation 'of  tL^T  8dVan?e  ?atio-  “ 
rotor  disk.  the  lnduced  velocity  over  the 

"SINC  BETAV1 

This  is  the  sine  of  the  local  flapping  angle. 

g.  "GOSC  BETA*) " 

This  is  the  cosine  of  the  local  flapping  angle. 

h.  "BETA  DOT" 

The  local  flapping  velocity  is  printed  in  degrees  per  second. 
2m  y.ar>iables  Which  Change  with  Radius 

a.  "RAD.  STA." 

rotoSrCra1dios1fo^eatch  lero^sm? ^  88  8  fractio"  the 

ea-h  aerodynamic  segment  of  the  blade. 

b.  "PHI" 

The  local  inflow  angle  is  printed  in  degrees  If  t-ho  «■ 
metric  incidence  is  zero,  this  is  the  angllof  attack.® 
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c.  "ALPHA” 

This  column  gives  the  local  angle  of  attack  for  each  blade  seg¬ 
ment. 

d.  "CL” 

The  local  lift  coefficients  are  printed. 

e.  "CD" 

The  local  drag  coefficients  are  printed. 

f.  ttMACH” 

The  local  Mach  number  is  given  for  each  blade  station. 

g.  "LOCAL  VI” 

The  local  values  of  induced  velocity  are  printed  in  feet  per 
second. 

h.  "LOC.  LAMBDA" 

The  symbol  lambda  is  usually  used  to  denote  the  inflow  ratio 
on  the  rotor.  Lambda,  then,  is  the  ratio  of  the  inflow  ve¬ 
locity  to  the  rotor  tip  speed.  For  a  stopped  rotor  lambda  is 
undefined,  so  the  output  in  this  column  is  the  inflow  velocity 
in  feet  per  second.  Positive  inflow  velocity  is  in  the  oppo¬ 
site  sense  of  positive  induced  velocity. 

i.  "UT" 

In  this  column  is  printed  the  velocity  component,  in  feet  per 
second,  which  is  perpendicular  to  the  radius  and  parallel  to 
the  local  zero  pitch  line  of  the  blade . 

j.  "UP" 

In  this  column  is  printed  the  Velocity  component,  in  feet  per 
second,  which  is  perpendicular  to  the  radius  and  perpendicular 
to  the  local  zero  pitch  line  of  the  blade. 


SECTION  IV 


DIAGNOSTIC  AMD  ERROR  MESSAGES 


All  of  the  messages  which  may  be  printed  out  from  the  computer 
program  and  which  are  considered  to  be  error  messages  are 
listed  below.  The  messages  are  in  alphabetical  order  with 
the  following  words  ignored:  M.R. ,  T.  R. ,  WING,  ELE,  FIN, 
COLLECTIVE  STICK,  F/A  CYCLIC  STICK,  LAT  CYCLIC  STICK,  PEDAL, 
MAIN,  TAIL,  WING,  ELEVATOR,  FIN/RUDDER. 

TVo  or  more  words  or  phrases  are  enclosed  in  brackets,  one 
above  the  other,  indicates  that  it  is  possible  to  have  either 
word  or  phrase,  but  only  one,  in  the  message  when  it  is 
printed  out.  An  underline  in  the  message  indicates  a  place 
for  a  numerical  value  in  the  message. 

After  each  message  is  the  name  of  the  subroutine  which 
printed  it  out.  The  next  statement  is  about  the  condition 
which  caused  the  message  to  be  printed.  Next  is  an  indica¬ 
tion  of  the  consequences  of  the  condition  followed  by  instruc¬ 
tions  to  the  user. 

1.  ALLEVIATION  DEVICE  FOR  ROTORS  BYPASSED  FOR  (^£1  ROTOR 

BECAUSE  ROTOR  DIAMETER  IS  TOO  SMALL  FOR  THIS  TIME 
INCREMENT  AND  VELOCITY. 

From  WAG 

The  analysis  in  WAG  assumes  a  minimum  number  of  data  points 
will  be  sampled  in  a  distance  traveled  which  is  calculated 
from  the  rotor  diameter.  This  message  indicates  that  the 
ratio,  VAt/ rotor  diameter,  is  too  large. 

WAG  is  bypassed  for  |^ail|  rotor* 


Make  At,  ZDEI.T1  or  ZDELT2  on  data  card  201,  smaller. 

2.  ALLEVIATION  DEVICE  FOR  WINGS  BYPASSED  BECAUSE  WING  CHORD 
IS  TOO  SMALL  FOR  THIS  TIME  INCREMENT  AND  VELOCITY. 

From  WAG 

The  analysis  contained  in  WAG  assumes  a  minimum  number  of 
data  points  will  be  sampled  in  a  distance  traveled  which  is 
calculated  from  the  wing  chord.  This  message  indicates  that 
the  ratio,  VAt/wing  chord,  is  too  large. 


WAG  is  bypassed  £or  wing. 

Make  At,  ZDELT1  or  ZDELT2  on  data  card  201,  smaller. 

3.  A  <_*_>  =  _ 

From  C0MS0L 

Matrix  of  complex  coefficients  of  linear  equations  in  the 
stability  analysis  is  singular. 

Program  continues. 

4.  CHECK  INPUT  FUSELAGE  INERTIAS.  THE  NUMBERS  INPUT  ARE 
PHYSICALLY  IMPOSSIBLE  AND  CANNOT  BE  HANDLED  BY  THIS 
PROGRAM. 

From  MNEM 

o 

This  message  indicates  that  Ix  Iz  -  1^  =  0. 

Problem  step  terminates. 

Change  the  input  data  for  Ix,  Iz  or  1^. 

5.  CHECK  PART  2  DATA  CARD  J  CODE  IS  _ 

From  IVAR 

A  value  for  J  on  data  card  211  has  been  input  for  which  an 
operation  is  not  defined. 

Problem  step  terminates. 

Change  the  card  indicated  by  the  message. 


M.R. 

6.  J  WING^  CLZ  FOR  I 
1  ELE  1 
FIN 


NORMAL 

REVERSED 


FLOW  HAS  BEEN  RESET  TO 


ALB  =  40  DEGREES. 


From  YFIX 


I?  J  in  the  Aerodynamic  Inputs  of  Section  II.  ALB 
-gure  11.  This  message  indicates  that  the  stall 


CLZ  is  YXX 
is  ag  in  Fi^ 

angle  for  Mach  number  =  0  is  too  large  (CLZ  is  too  large)  and 
thus  has  been  reduced. 


7.  DATA  ERROR  ..  NPART  = 


From  C81 

The  control  program,  C81  read  an  Illegal  value  of  NPART  on 
data  card  01.  This  error  most  commonly  occurs  after  another 
error  has  interrupted  the  normal  sequence  of  events  by  ter¬ 
minating  the  problem  step. 

Program  execution  terminates. 


\  wing} 

I  ELE  ( 
\FIN  / 


DRAG  DIVERGENCE  MACH  NUMBER  HAS  BEEN  RESET 


TO  _ . 

From  YFIX 

YXX(l)  was  input  greater  than  or  equal  to  1.  This  is  a 
warning  message. 

9.  EXCESSIVE  ANGLE  OF  ATTACK  FOR  N  = 


From  CLCD 


Subroutine  CLCD  was  entered  with  the  angle  of  attack  of  the 

(main  rotor\ 
tail  rotor j 

wing  >  greater  than  20  radians, 
elevator  1 
f in/rudder/ 

Problem  step  terminates. 

10.  FLAPPING  ANGLES  CANNOT  BE  ESTIMATED  FOR  A  STOPPED  ROTOR. 
From  MNEM 

The  input  data  are  contradictory.  XFC  (22)  =0  and  either 
XFC  (25)  =  0  or  XMR  (13)  =  0  or  both. 


Problem  step  terminates. 


Change  the  input  data  to  either  XFC(22)  =  0  or  XFC(25)  =  0 
and  XMR( 13)  =  0. 

11.  FUNCTION  VALUE  IS  ZERO  FOR  ALL  VALUES  OF  X. 

From  R00A. 

Hie  determinant  of  one  of  the  characteristic  coefficient 
matrices  in  the  stability  analysis  is  identically  zero. 

Program  continues. 

Check  coefficients  for  realistic  values. 

12.  INDUCED  VELOCITY  SET  TO  0.  CALCULATIONS  NON  CONVERGENT. 
From  VI ND 

The  iterative  loop  which  calculates  induced  velocity  for  a 
constant  velocity  and  rotor  thrust  has  not  converged  in  100 
iterations. 

Warning  message  only. 

13.  IPSN  INDICATED  NOT  ON  LIBRARY. 

From  C81L 

In  an  operation  with  NPART  =  8,  NSCALE  =  2  on  Card  01,  the 
IPSN  input  on  card  02  does  not  match  any  IPSN  on  the  file 
tape. 

Problem  step  terminates. 

Check  input  IPSN  and  list  of  IPSN' s  on  the  file  tape. 


14. 


MACH  NUMBER  FOR  LOWER  BOUND  OF  SUPERSONIC  REGION 


HAS  BEEN  RESET  TO  _ . 

From  YFIX 

YXX(2)  was  input  less  than  or  equal  to  1. 
This  is  a  warning  message. 
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15.  NEXT  ROOT  GREATER  THAN  _  RADIANS  INCOMPLETE 

FUNCTION,  RESIDUE  F(S)  TJEJKTNING. 

From  ROOA 

One  of  the  characteristic  coefficient  matrices  in  the 
stability  analysis  is  ill-conditioned. 

Program  continues. 

Check  coefficients  for  realistic  values. 

16.  NSCALE  HAS  ILLEGAL  VALUE. 

From  C81L 

On  data  card  01  with  NPART  =  8,  NSCALE  has  a  value  for  which 
no  operation  is  defined. 

Problem  step  terminates. 

17.  THE  PARTIAI  DERIVATIVE  MATRIX  IS  SINGULAR.  THIS  IS 
PROBABLY  CAUSED  BY  ONE  OF  1HE  CONTROLS  BEING  UNCONNECTED. 

From  I1RIM 

During  the  TRIM  procedure,  a  singular  partial  derivative 
matrix  occurred.  The  usual  cause  is  an  error  in  the  input 
data  for  one  of  the  controls.  Previous  matrices,  if  ahy, 
should  be  examined  for  a  near-zero  row  or  column  to  help 
locate  the  cause. 

Problem  step  terminates. 

18.  THE  PHASE  ANGLE  DIFFERENCE  BETWEEN  AND 

IS  A  MULTIPLE  OF  180  DEGREES.  THEREFORE  tfO  VARIXBl'E  'CAM  ' 
BE  EXPRESSED  AS  A  LINEAR  FUNCTION  OF  THEM. 

From  CURVET 

The  Vector  analysis  section  of  the  program  where  the  coef¬ 
ficients  in  the  expression  A  =  K  B  *  B  +  KC  *  C  +  D  are 
derived  has  failed  because  of  the  linear  dependency  of  B  and 
C. 

Program  goes  to  next  set  of  variables. 
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! COLLECTIVE  STICK  \ 

lat  cyclic  stick!  position  exceeds  stops 

PEDAL  / 

PERCENT  FULL  THROW  COMPUTED. 


From  TRIM 


The  program  found  a  stable  flight  condition,  but  one  of  the 
control  positions  necessary  is  outside  the  allowable  range. 

Change  flight  conditions,  configuration,  or  allowable  control 
settings. 

20.  POWER  REQUIRED  FOR  TRIM  CONDITION  EXCEEDS  POWER  AVAILABLE. 
From  PARA 


The  power  required  to  maintain  constant  rpm  at  the  trim 
condition  is  greater  than  the  engine  is  capable  of  producing. 


Check  for  rotor  stall,  practicality  of  flight  conditions,  or 
increase  engine  power ? 


21. 


MAIN 

TAIL 


ROTOR  BALANCE  EXCEEDS  ALLOWABLE  ERROR. 


From  I TOOT 


The  iteration  loop  in  the  rotor  analysis  which  balances  the 
rotor  flapping  moments  was  activated  but  could  not  balance 
the  rotor  in  the  number  of  iterations  allowed. 


Problem  step  terminates. 


Check  configuration,  flight  regime,  and  spatial  orientation 
for  compatibility. 


22. 


(MAIN 

(TAIL 


ROTOR  FLAP  CORRECTION  EXCEEDS  90  DEGREES 


From  ITROT 


The  iteration  loop  in  the  rotor  analysis  which  balances  the 
rotor  flapping  moments  was  activated  and  computed  a  cor¬ 
rection  to  a  flapping  angle  which  exceeded  90  degrees. 


Problem  step  terminates. 

Check  configuration,  flight  regime,  and  spatial  orientation 
for  compatibility. 


23.  SHIP  CONTACTS  GROUND 
From  VIND 

Altitude  has  become  negative. 

Problem  stop  terminates. 

Find  out  why  nhip  lost  altitude  and  correct. 

24.  SOLUTION  EXCEEDS  MAXIMUM  NUMBER  OF  ROOTS  INCOMPLETE 
FUNCTION,  RESIDUE  F(S)  REMAINING. 

From  ROOA 

One  of  the  characteristic  coefficient  matrices  in  the 
stability  analysis  is  ill-conditioned. 

Program  continues. 

Check  coefficients  for  realistic  values. 

(  WING  ) 

25.  {ELEVATOR  !  STALLED  AT  DEGREES. 

( FIN/RUDDER) 

CL  =  _  CD  =  _ 

From  CLCD 

This  is  a  warning  that  the  surface  is  stalling  out. 

26.  WARNING,  THE  PARTIAL  DERIVATIVE  MATRIX  MAY  BE  IN  ERROR. 
From  ITROT 

In  the  rotor  analysis  the  iteration  loop  which  balances  the 
rotor  flapping  moments  and  the  thrust- induced  velocity 
iteration  loop  are  both  activated.  While  each  is  able  to 
converge  separately,  they  have  not  been  able  to  converge 
together. 

Warning  message. 

Exercise  care  in  use  of  the  partial  derivative  matrix  follow¬ 
ing  this  message. 
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Figure  2.  Parameter  Sweep  With  Trim 


Figure  6.  Maneuver  With  Stability  Analyst 
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Figure  11.  Sample  Data  Set  for  Single  Rotor  Helicopter  Case 
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M103MI  SAMPLE  DATA  PM  HORIZONTAL  STOP  ANO  POLC  HELICOPTER 

UIN6  AREA  *  500. 0  SO.  ET.  NINO  CM  -  0.0005 

CHECK  OECK  EUR  HORRI ZOMTAL  EOLO  -RUOOER  LINKED  TO  PEDAL 
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Figure  12*  Sample  Data  Set  for  Horizontal  Stop 
and  Fold  Helicopter  Case 
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Figure  13.  Sample  Data  Set  for  Tilting  Prop-Rotor 
In  Helicopter  Mode 
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Figure  13,  Concluded 


O  O  OOO  0000*00  OOO  0000*00  ooo  o 

80  OOO  0000*00  OOO  0000*00  ooo  o 

O  OOO  0000*00  OOO  0000*00  ooo  P 

o  O  OOO  *000*00  OOO  *000*00  ooo  O 

s.  S  888  *:0.88*88  °.88  T?°.8*.§8  §28  8 

30*00  O  •  »OK^**«  *000  O  ••©>»♦**♦  *000  OM  i  *00  00*00 

•  *4  •  •  N«4N  •  i«*  •  •  •  N«N  *AN  N  •  •  **  •  •  •  •  *•*«  •  •  •  •  *4  •  • 

O  OO  O  OO  OOO  O  OO  OOO  OO  OO  OO  OO 


0000*00  OOO  0000*00 
0000*00  OOO  OQOO*Op 
0000*00  OOO  0000*00 
*000*00  OOO  *000*00 

°*88  *.*.*.%*$ 8 


* 

* 

ooo  •  o 

(  •  1*4  1 

ooo  o 


5  888  888888 

*  "88  S88SS8 

888  'O'OSSSS 


OOO  •  tO«4N 

O  •  *oo*»  •*  «  N  O  O  O 
O  OOOOOO 


888  §88888 
°.§8  §88888 
888  .*823° 


O  •  *OOK  **#NOOO 

•*»*N  •  . 

O  OOOOOO 


0«*40  OO 
•  ••••• 
OOOO OO 


»ooo  • 

I  *  *  *  * 

ooo 


88  §888888  S  88  §888888  S 

8§  S28|§88  §  82  SS|8888  2 

O  4  NK05o*0  *  04  Nh»  5  O  Oir  O  * 

•  O  •#  *  *  *  *  •  O  *  *4  *  *  *  *  *  •  o  * 


1  A 
8§8  § 

ill  i 


•  4  *  •  •**«•  o  * 

*  •  O  *  *  *4  *44  AOKO 
«4  4  I*4«HM4  ••  •  •  •  • 

o  OOOOOO 


•  4  *••**«  o  * 

4  •  O  •  *  •  4  4  *  O  h*  O 

*4  4  *>4MM  4  •  •  •  ••• 
O  OOOOOO 


•  *  *  O KO 
<#••••• 
OOO  OO 


?  3 


2  388858  §884?  8888888 

3  §88853  5883  8S888S8 

mo  •**«  <*o  •  o«40*ooo 

^  *•*»•*•  •  •  *  •  •  *N^4  O 

3  34*4440  2NO0O*4MOKmMK 


§888  §88888! 


34*4440  ZN0(*0*4N0K^MN0  0  JN040- 

lk«4*4  •  *4  M  *4  *  *4  •  44*4  «  •  •  f  •  •  ***44*4  • 

I  O  4  O  OOOOOO  <  O 


ooSS  sssssss 

«*©  •  ONOfOOO 
•  *4  •  •  •  4*  rn*  O 

JN040*4N0Kr  4*-0  0 
•-***•4  *  *  *  *4  •••••• 

4  O  OOOOOO 


3  - 

5  I 
28  8 
*8  8 
°.  s 

i0.- 

oo 

I 


8  §8  888§8o888§§ 

•  *  o  oo*Oooooo** 


888888888 §§ 
00*000000** 
00*0*0  0  O  o  *  * 


O  *0  00*0*0000**  00*0*0000** 

8  §4  85i,S8.‘S*.*.|sift  'Hwmztz 

0410. 4*  404444*  *  *  *  *OQ  404*444  **4  *00 

•  *  •  •  *4  *4  *4  «4  *4  4  A  *4  >4  •  •  *4  •  •  *4*4*4*44A*4*4  •  •  *4  •  • 

O  OO  I  oo  oo  •  oo  oo 


o 

A 

88585 
00*0* 
O  •  *  o  * 

J5S21S 

? 


*0044^  *4A  •  O  W 
t  •  *4A44  »  i*4  i  « 

OO  O  O  O  O 


O  >004  44  •  *  *  *041  *  •  O  •  O  i 

4  *4  •  *44*44  •  *4  •  •  44  •  *4  • 

900  oooo  o  o  < 


8858 


•  0*0 
00*4 
O*  •* 
OO*  •  < 


88  8838885588  88  8888§8$888 

88  8S88585S88  §8  8288S35S88 


OO  O  •  •  •*i*»4*oO 
•  OO  *440  1440*4 
**  •  *  4*4*  •  •  •  •  • 

o  o  o  o  o  o 


OO  O  •  •  (44*4  Odu 
•  OO  *440  *440*4 
44  *44*4*  •  •  •  •  • 

o  o  o  o  o  o 


118 


. 

_  _ 


o  oo  oo 
•  •  •  •  • 
o  oooo 


o 

0  0  0-40 
•  •  •  •  • 
ooooo 


S  §8! 
§  58! 

3»  •  « 

IW4*  I 


O  O  •  J*4 

•  •  41  Uf  « 

o  o 


S  5 
A  A 

§558  I 
0**0  * 
0**0  * 

3J8?  1 


88  I 

oo  * 


oo  * 

°.s  s 

oo  .o» 


88  8885588  8  §8  8885588  |  85  8  5 

88  8285888  §  8§  §285538  8  85  8  5 

O  •  *40*4004  O  O  •  *KO*«*ON  O  •  4  O  * 

•  *  •  •  •4*04  O  •*  •  •  *4*04  O  *0  O  * 


§88  § 
§88  8 
^88  8 


88 

88 

ii 

888 

88 

§8 

S§§ 

§5 

•  o 
o 
►-  m 

ii 

<4*00 

OOO 
•  UN 

0*0 

ooo  oo 


_ ■ 


Figure  14.  Sample  Data  Set  for  Tilting  Prop-Rotor 
In  Airplane  Mode 
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Precone  Angle 
Degrees 

Rotor  Torque 
Foot -Pounds 

Rotor  Speed 
RPM 

Rotor  Thrust 
Pounds 


Maneuver  Time  Seconds 


Figure  16.  Time  Histories  for  Tilt-Forward 

Trail-Aft  Case 


Figure  18.  Input  Data  for 


Figure  20.  Final  Trim  Page 
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Figure  22.  Typical  Maneuver  Page 


Figure  26.  Stability  Partial  Derivative  Matrices 
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Figure  29.  Harmonic  Analysis  Results  for 
Low  Frequency  Case 
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Figure  30.  Amplitude  and  Phase  Angle  Comparison  for 
Vector  Analysis  at  Low  Frequency 


COMP  It  CO  ll/tt/M 
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Figure  31.  One  Variable  as  a  Function  of  Two  Other 
Variables  for  Low  Frequency  Case 
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Rotor  Airload  Data  Output 


|  TABLE  I.  PLOT  INDEX  1 

Group 

Number 

Description 

FIRST 

I 

Q  VELOCITY,  TPP1,  DEG/SEC 

ROTOR 

2 

P  VELOCITY,  IPP1,  DEG/SEC 

GROUP 

3 

U  VELOCITY,  MAST1  AXES,  FT/SEC 

4 

OMEGA-DOT,  R0T0R1,  DEG/SEC/SEC 

5 

BETA-  2DOT,  BLADE  1 ,  ROTOR1 , DEG/SEC/SBC 

6 

R0T0R1  THRUST,  LB 

7 

F/A  FLAPPING,  MASTl/TPPl,  DEG 

8 

LATERAL  FLAPPING,  MASTl/TPPl,  DEG' 

9 

V  VELOCITY,  MAST1  AXES,  FT/SEC 

10 

OMEGA,  ROTOR1, DEG/SEC 

II 

BETA-DOT,  BLADE1,  ROTOR1,  DEG/SEC 

12 

ROTOR1  H- FORCE,  LB 

13 

W  VELOCITY,  MAST!  AXES,  FT/SEC 

14 

AZIMUTH  LOC. ,  BLADE1 ,  R0T0R1,  DEG 

15 

BETA,  BLADE1,  ROTOR1,  DBG 

16 

ROTOR1  Y- FORCE,  LB 

17 

ROTOR1  COLLEC.  FROM  CONTROLS,  DEG 

18 

ROTOR1  F/A  CYC.  FROM  CONTROLS,  DEG 

19 

ROTOR1  LAT  CYC.  FROM  CONTROLS,  DEG 

20 

ROTOR1  CONING,  DEG 

21 

MAST1  TILT,  DEG 

22 

ROTOR1  TORQUE,  FT- LB 

23 

ROTOR1  OTHER  COLLEC. ,  DEG 

24 

ROTOR1  OTHER  F/A  CYC.  ,  DEG 

25 

ROTOR1  OTHER  LAT  CYC. ,  DEG 

26 

ROTOR1  INDUCED  VELOCITY,  FT/SEC 

27 

ROTOR1  RPM 

28 

R0T0R1  HORSEPOWER 

29 

ROTOR1  UPPER  FLAPPING  LIMIT,  DEG 

30 

ROTOR1  F/A  HUBSPRING,  FT-LB/DEG 

31 

ROTOR1  TOTAL  COLLECTIVE,  DEG 

32 

ROTOR1  TOTAL  F/A  CYCLIC,  DEG 

33 

ROTOR1  TOTAL  LAT  CYCLIC,  DEG 

34 

ROTOR1  LOWER  FLAPPING  LIMIT,  DEG 

35 

ROTOR1  LAT  HUBSPRING,  FT-LB/DEG 

SECOND 

36 

Q  VELOCITY,  TPP2,  DEG/SEC 

ROTOR 

37 

P  VELOCITY,  TPP2 ,  DEG/SEC 

GROUP 

38 

U  VELOCITY,  MAST2  AXES,  FT/SEC 

39 

OMEGA-DOT,  R0T0R2,  DEG/SEC/SEC 

40 

BETA- 2DOT, BLADE1  ,ROTOR2 , DEG/SEC/SEC 

41 

ROTOR2  THRUST,  LB 

42 

F/A  FLAPPING,  MAST2/TPP2 ,  DEG 

43 

LATERAL  FLAPPING,  MAST2/TPP2,  DEG 

44 

V  VELOCITY,  MAST2  AXES,  FT/SEC 

45 

OMEGA,  RDTOR2,  DEG/SEC 

46 

BETA-DOT.  BLADE1.  ROTOR2.  DEG/SEC 
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TABLE  I.  Continued 

Group 

Number 

Description 

SECOND 

47 

ROTOR2  H- FORCE,  LB 

ROTOR 

48 

W  VELOCITY,  MAST2  AXES,  FT/SEC 

GROUP 

49 

AZIMUTH  LOC. ,  BLADE!,  ROTOR2,  DEG 

(Cont'd) 

50 

BETA,  BLADE!,  ROTOR 2,  DEG 

51 

R0TOR2  Y- FORCE,  LB 

52 

R0TOR2  COLLEC.  FROM  CONTROLS,  DEG 

53 

ROTOR2  F/A  CYC.  FROM  CONTROLS,  DEG 

54 

R0TOR2  LAT  CYC.  FROM  CONTROLS,  DEG 

55 

R0TOR2  CONING,  DBG 

56 

MAST2  TILT,  DEG 

57 

ROTOR2  TORQUE,  FT- LB 

58 

ROTOR2  OTHER  COLLEC. ,  DEG 

59 

R0TOR2  OTHER  F/A  CYC. ,  DEG 

60 

R0TOR2  OTHER  LAT  CYC. ,  DEG 

61 

ROTOR2  INDUCED  VELOCITY,  FT/SEC 

62 

R0T0R2  RPM 

63 

ROTOR2  HORSEPOWER 

64 

R0TOR2  UPPER  FLAPPING  LIMIT,  DEG 

65 

R0T0R2  F/A  HUBSPIRNG,  FT-LB/DEG 

66 

R0TOR2  TOTAL  COLLECTIVE,  DEG 

67 

R0TOR2  TOTAL  F/A  CYCLIC,  DEG 

68 

R0TOR2  TOTAL  LAT  CYCLIC,  DEG 

69 

R0T0R2  LOWER  FLAPPING  LIMIT,  DEG 

70 

ROTOR 2  LAT  HUBSPRING,  FT-LB/DEG 

FIXED 

71 

X-COMP  VELOCITY,  FIXED  AXES,  FT/SEC 

REFERENCE 

72 

Y-CCMP  VELOCITY,  FIXED  AXES,  FT/SEC 

GROUP 

73 

Z-COMP  VELOCITY,  FIXED  AXES,  FT/SEC 

74 

TOTAL  DISTANCE  FLOWN,  FT 

75 

AIR  SPEED,  KTS 

76 

HEADING  ANGLE,  DEG 

77 

X-COMP  DISP. ,  FIXED  AXES,  FT 

78 

Y-COMP  DISP.,  FIXED  AXES,  FT 

79 

Z-COMP  DISP. ,  FIXED  AXES,  FT 

80 

ALTITUDE,  FT 

81 

GROUND  SPEED,  KTS 

82 

CLIMB  ANGLE,  DEG 

FUSELAGE 

83 

U-DOT  ACCEL. ,  BODY  AXES,  FT/SEC/SEC 

GROUP 

84 

V-DOT  ACCEL.,  BODY  AXES,  FT/SEC/SEC 

85 

W-DOT  ACCEL. ,  BODY  AXES,  FT/SEC/SEC 

86 

P-DOT  ACCEL,  BODY  AXES,  DBG/SEC/SEC 

87 

Q-DOT  ACCEL,  BODY  AXES,  DEG/SEC/SEC 

88 

R-DOT  ACCEL,  BODY  AXES,  DEG/SEC/SEC 

89 

COLLEC.  BOBWT.  ACCEL.,  DEG/SEC/SEC 

90 

U  VELOCITY,  BODY  AXES,  FT/SEC 

91 

V  VELOCITY.  BODY  AXES.  FT/SEC 

. 


Wt 


Group  Humber 


TABLE  I.  Continued 
Description 


FUSELAGE 

92 

W  VELOCITY,  BODY  AXES,  FT/SEC 

GROUP 

93 

P  VELOCITY,  BODY  AXES,  DEG/SEC 

(Cont’d) 

94 

Q  VELOCITY,  BODY  AXES,  DEG/SBC 

95 

R  VELOCITY,  BODY  AXES,  DEG/SEC 

96 

COLLEC.  BOBWT.  VELOCITY,  DEG/SEC 

97 

YAW  VELOCITY,  FIXED/BODY,  DEG/SEC 

98 

PITCH  VELOCITY,  FIXED/BODY,  DEC/SEC 

99 

ROLL  VELOCITY,  FIXED/BODY,  DEG/SEC 

100 

YAW  ANGLE,  FIXED/BODY,  DBG 

101 

PITCH  ANGLE,  FIXED/BODY,  DEG 

102 

ROLL  ANGLE,  FIXED/BODY,  DEG 

103 

COLLECTIVE  STICK  POSITION,  PCT 

104 

F/A  CYCLIC  STICK  POSITION,  PCT 

105 

LEFT  WING  ANGLE  OF  ATTACK,  DEG 

106 

RIGHT  WING  ANGLE  OF  ATTACK,  DEG 

107 

ELEVATOR  ANGLE  OF  ATTACK,  DEG 

108 

FIN  ANGLE  OF  ATTACK,  DEG 

109 

BODY  YAW  WRT.  FLIGHT  PATH,  DEG 

110 

C.G.  STATION  LINE  LOCATION,  IN. 

111 

X-COMP  GUST  VEL. ,  BODY  AXES,  FT/SEC 

112 

X-COMP  G-S,  BODY  AXES 

113 

LATERAL  CYCLIC  STICK  POSITION,  PCT 

114 

LEFT  WING  COEFFICIENT  OF  LIFT 

115 

RIGHT  WING  COEFFICIENT  OF  LIFT 

116 

ELEVATOR  COEFFICIENT  OF  LIFT 

117 

FIN  COEFFICIENT  OF  LIFT 

118 

BODY  PITCH  WRT.  FLIGHT  PATH,  DEG 

119 

C.G.  BUTT  LINE  LOCATION,  IN. 

120 

Y-CGMP  GUST  VEL.,  BODY  AXES,  FT/SEC 

121 

Y-COHP  G-S,  BODY  AXES 

122 

PEDAL  POSITION,  PCT 

123 

LEFT  WING  COEFFICIENT  OF  DRAG 

124 

RIGHT  WING  COEFFICIENT  OF  DRAG 

125 

ELEVATOR  COEFFICIENT  OF  DRAG 

126 

FIN  COEFFICIENT  OF  DRAG 

127 

C.G.  WATER  LINE  LOCATION,  IN. 

128 

Z-COMP  GUST  VEL.,  BODY  AXES,  FT/SBC 

129 

Z-COMP  G-S,  BODY  AXES 

130 

RIGHT/CENTER  JET  THRUST,  LB 

131 

ENGINE  TORQUE  SUPPLIED,  FT-LB 

132 

TOTAL  HORSEPOWER  REQUIRED 

133 

LEFT  JET  THRUST,  LB 

134 

SHAFT  HORSEPOWER 

135 

ROTOR  BRAKE  TORQUE  APPLIED,  FT-LB 

Group  Number 


TABLE  I.  Continued 
Description 


AZIMUTH 

136 

AZIMUTH  LOC. ,  BLADE1,  R0T0R1,  DEG 

LOCATION 

137 

AZIMUTH  LOC. ,  BLADE 2,  R0T0R1 ,  DEG 

AND 

138 

AZIMUTH  LOC. ,  BLADE 3,  ROTORl,  DEG 

FLAPPING 

139 

AZIMUTH  LOC. ,  BLADE 4,  ROTOR! ,  DEG 

GROUP 

140 

AZIMUTH  LOC. ,  BLADE5,  R0T0R1,  DEG 

141 

AZIMUTH  LOC. ,  BLADE6,  ROTORL ,  DBG 

142 

AZIMUTH  LOC. ,  BLADE 7,  ROTORL,  DEG 

143 

BETA-  2D0T,  BLADEL ,  ROTORL ,  DEG/SEC/SEC 

144 

BETA-  2  DOT,  BLADE2 ,  R0TCR1 ,  DEG/SBC/SEC 

145 

BETA-  2D0T,  BLADE3 ,  R0T0R1 ,  DEG/SEC/SEC 

146 

BETA-  2D0T ,  BLADE  4,  ROTORL ,  DEG/SEC/SEC 

147 

BETA-  2D0T,  BLADE5 ,  ROTORl ,  DEG/SBC/SBC 

148 

BETA-  2D0T,  BLADE6 ,  RO  TORI ,  DBG/SEC/SEC 

149 

BETA-  2D0T ,  BLADE7 ,  ROTORl ,  DEG/S  EC/SEC 

150 

BETA- DOT,  B  LAD  El ,  ROTORl,  DEG/SBC 

151 

BETA-DOT,  BLADE2,  ROTORl,  DEG/SEC 

152 

BETA-DOT,  BLADE3,  ROTORl,  DEG/SEC 

153 

BETA-DOT,  BLADE4,  ROTORl,  DEG/SEc 

154 

BETA-DOT,  BLADE5,  ROTORl,  DEG/SEC 

155 

BETA-DOT,  BLADE6,  ROTORl,  DEG/SEC 

156 

BETA-DOT,  BLADE7,  ROTORl,  DEG/SEc 

157 

BETA,  BLADE1,  ROTORl,  DEG 

158 

BETA,  BIADE2,  ROTORl,  DEG 

159 

BETA,  BLADE3,  ROTORl,  DEG 

160 

BETA,  BLADE4,  ROTORl,  DEG 

161 

BETA,  BLADE5,  ROTORl,  DEG 

162 

BETA,  BLADE6,  ROTORl,  DEG 

163 

BETA,  BLADE7,  ROTORl,  DEG 

164 

AZIMUTH  LOC.,  BLADE1,  ROTOR2,  DEG 

165 

AZIMUTH  LOC.,  BLADE2,  ROTOR2,  DEG 

166 

AZIMUTH  LOC.,  BLADE3,  ROTOR 2 ,  DEG 

167 

AZIMUTH  LOC.,  BLADE4,  R0T0R2,  DEG 

168 

AZIMUTH  LOC.,  BLADE5,  ROTOR 2 ,  DEG 

169 

AZIMUTH  LOC.,  BLADE6,  ROTOR 2,  DEG 

170 

AZIMUTH  LOC.,  BLADE7,  ROTOR2,  DEG 

171 

BETA-  2D0T,  BLADE1 ,  R0TOR2 ,  DEG/SEC/SEC 

172 

BETA-  2D0T ,  BLADE2 ,  RO  T0R2 ,  DEG/S  EC/SEC 

173 

BETA-  2D0T,  BLADE3 ,  R0T0R2 ,  DEG/  SEC/SEC 

174 

BETA-  2D0T,  BLADE4 ,  ROTOR2 ,  DEG/SEC/SEC 

175 

BETA-  2  DOT,  BLADE5 ,  R0T0R2 ,  DEG/S  EC/SEC 

176 

BETA-  2  DOT,  BLADE6  ,ROTOR2 ,  DEG/SBC/SEC 

177 

BETA- 2DOT,BLADE7,ROTOR2, DEG/SEC/SEC  ; 

178 

BETA-DOT,  BLADEL,  R070R2,  DEG/SEC 

179 

BETA-DOT,  BLADE2,  R0T0R2,  DEG/SEC 

180 

BETA- DOT,  BLADE 3,  R0T0R2,  DEG/SEC 

181 

BETA-DOT,  BLADE4,  R0T0R2,  DEG/SEC 

182 

BETA- DOT.  BLADE5,  ROTOR2,  DEG/SEC 

TABLE  1.  Continued 

Group 

Number 

Description 

AZIMUTH 

183 

BETA-DOT,  BLADE6,  R0T0R2,  DEG/SEC 

LOCATION 

184 

BETA- DOT,  BLADE7,  ROTOR 2,  DEG/SEC 

AND 

185 

BETA,  BLADE1,  R0T0R2,  DEG 

FLAPPING 

186 

BETA,  BLADE2,  R0T0R2,  DEG 

GROUP 

187 

BETA,  BLADE3,  ROTOR2,  DEG 

(Cont’d) 

188 

BETA,  BLADE4,  R0T0R2,  DEG 

189 

BETA,  BLADE5,  R0T0R2,  DEG 

190 

BETA,  BLADE6,  R0T0R2,  DEG 

191 

BETA,  BLADE7,  R0T0R2,  DEG 

FORCE  AND 

192 

TOTAL  X-FORCE  ON  C.G. ,  LB 

MOMENT 

193 

X- FORCE  FROM  RIGHT  WING,  LB 

GROUP 

194 

X-FORCE  FROM  LEFT  WING,  LB 

195 

X-FORCE  FROM  ELEVATOR,  LB 

196 

X-FORCE  FROM  FUSELAGE,  LB 

197 

X-FORCE  FROM  RIGHT  JET,  LB 

198 

X-FORCE  FROM  LEFT/CENTER  JET,  LB 

199 

X-FORCE  FROM  R0T0R1,  LB 

200 

X-FORCE  FROM  R070R2,  LB 

201 

X-FORCE  FROM  WEAPON  FIRE,  LB 

202 

X-FORCE  FROM  FIN,  LB 

203 

X-FORCE  FROM  WEIGHT,  LB 

204 

TOTAL  Y-FORCE  ON  C.G. ,  LB 

205 

Y-FORCE  FROM  FUSELAGE,  LB 

206 

Y-FORCE  FROM  RIGHT  JET,  LB 

207 

Y-FORCE  FROM  LEFT/CENTER  JET,  LB 

208 

Y-FORCE  FROM  R0T0R1,  LB 

209 

Y-FORCE  FROM  R0T0R2,  LB 

210 

Y-FORCE  FROM  WEAPON  FIRE,  LB 

211 

Y-FORCE  FROM  FIN,  LB 

212 

Y-FORCE  FROM  WEIGHT,  LB 

213 

TOTAL  Z-FORCE  ON  C.G. ,  LB 

214 

Z -FORCE  FROM  RIGHT  WING,  LB 

215 

Z-FORCE  FROM  LEFT  WING,  LB 

216 

Z-FORCE  FROM  ELEVATOR,  LB 

217 

Z-FORCE  FROM  FUSELAGE,  LB 

218 

Z-FORCE  FROM  RIGHT  JET,  LB 

219 

Z-FORCE  FROM  LEFT/C  ENTER  JET,  LB 

220 

Z-FORCE  FROM  R0T0R1,  LB 

221 

Z-FORCE  FROM  ROTOR2,  LB 

222 

Z-FORCE  FROM  WEAPON  FIRE,  LB 

223 

Z-FORCE  FROM  WEIGHT,  LB 

224 

TOTAL  ROLL  MOM  ON  C.G. ,  FT-LB 

225 

ROLL  MOM  FROM  RIGHT  WING,  LB 

226 

ROLL  MOM  FROM  LEFT  WING,  FT-LB 

227 

ROLL  MOM  FROM  ELEVATOR,  FT-LB 

228 

ROLL  MOM  FROM  FUSELAGE.  FT-LB 
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Group  Number 


TABLE  I.  Continued 
Description 


FORCE  AND 

229 

ROLL  MOM  FROM  RIGHT  JET,  FT- LB 

MOMENT 

230 

ROLL  MOM  FROM  LEFT/C ENTR  JET,  FT-LB 

GROUP 

231 

ROLL  MOM  FROM  ROTOR1  FORCES,  FT-LB 

(Cont'd) 

232 

ROLL  MOM  FROM  R0TOR2  FORCES,  FT-LB 

233 

ROLL  MOM  FROM  WEAPON  FIRE,  FT-LB 

234 

ROLL  MOM  FROM  FIN,  FT-LB 

235 

ROLL  MOM  FROM  ROTOR!  TORQUE,  FT-LB 

236 

ROLL  MOM  FROM  R0TOR2  TORQUE,  FT-LB 

237 

TOTAL  PITCH  MOM  ON  C.G. ,  FT-LB 

238 

PITCH  MOM  FROM  RIGHT  WING,  FT-LB 

239 

PITCH  MOM  FROM  LEFT  WING,  FT-LB 

240 

PITCH  MOM  FROM  ELEVATOR,  FT-LB 

241 

PITCH  MOM  FROM  FUSELAGE,  FT-LB 

242 

PITCH  MOM  FROM  RIGHT  JET,  FT-LB 

243 

PITCH  MOM  FROM  LEFT/CENT  JET,  FT-LB 

244 

PITCH  MOM  FROM  R0T0R1  FORCES,  FT-LB 

245 

PITCH  MOM  FROM  ROTOR2  FORCES,  FT-LB  | 

246 

PITCH  MOM  FROM  WEAPON  FIRE,  FT-LB 

247 

PITCH  MOM  FROM  FIN,  FT-LB 

248 

PITCH  .MOM  FROM  ROTOR1  TORQUE,  FT-LB 

249 

PITCH  MOM  FROM  RDT0R2  TORQUE,  FT-LB 

250 

TOTAL  YAW  MOM  ON  C.G.,  FT-LB 

251 

YAW  MOM  FROM  RIGHT  WING,  FT-LB 

252 

YAW  MOM  FROM  LEFT  WING,  FT-LB 

253 

YAW  MOM  FROM  ELEVATOR,  FT-LB 

254 

YAW  MOM  FROM  FUSELAGE,  FT-LB 

255 

YAW  MOM  FROM  RIGHT  JET,  FT-LB 

256 

YAW  MOM  FROM  LEFT/CENT  JET,  FT-LB 

257 

YAW  MOM  FROM  ROTORl  FORCES,  FT-LB 

258 

YAW  MOM  FROM  ROTOR2  FORCES,  FT-LB 

259 

YAW  MOM  FROM  WEAPON  FIRE,  FT-LB 

260 

YAW  MOM  FROM  FIN,  FT-LB 

261 

YAW  MOM  FROM  ROTORl  TORQUE,  FT-LB 

262 

YAW  MOM  FROM  ROTOR 2  TORQUE,  FT-LB 

263 

ROTORl  F/A  FLAPPING  MOMENT,  TT-LB 

264 

ROTORl  LAT  FLAPPING  MOMENT,  FT-LB 

265 

ROTOR 2  F/A  FLAPPING  MOMENT,  FT-LB 

266 

ROTOR2  LAT  FLAPPING  MOMENT,  FT-LB 
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TABLE  II.  CONTROL  TO  SURFACES 


Configuration 

Control 

Surface 

Single  rotor 

Collective  stick 

Main  rotor  collec¬ 
tive 

Wing  incidence 
Elevator  incidence 
Jet  thrust 

F/A  cyclic  stick 

Main  rotor  F/A 
cyclic 

Elevator  incidence 

Lateral  cyclic  stick 

Main  rotor  lateral 
cyclic 

Aileron 

Pedal 

Tail  rotor  collec¬ 
tive 

Rudder  incidence 

Tandem 

Collective  stick 

Fwd  rotor  collec¬ 
tive 

Aft  rotor  collec¬ 
tive 

Wing  incidence 
Elevator  Incidence 
Jet  thrust 

F/A  cyclic  stick 

Fwd  rotor  F/A  cyclic 
Aft  rotor  F/A  cyclic 
Fwd  rotor  collective 
Aft  rotor  collective 
Elevator  incidence 

Lateral  cyclic  stick 

Fwd  rotor  lateral 
cyclic 

Aft  rotor  lateral 
cyclic 

Aileron 

Pedal 

Rudder  incidence 

Fwd  rotor  lateral 
cyclic 

Aft  rotor  lateral 
cyclic 

TABLE  II. Continued 


Configuration 

Control 

Surface 

Composite 

Collective  stick 

Rt  rotor  collective. 
Left  rotor  collec¬ 
tive 

Wing  incidence 
Elevator  incidence 

Jet  thrust 

F/A  cyclic  stick 

Rt  rotor  F/A  cyclic 
Left  rotor  F/A  cyclic 
Rt  rotor  lateral 
cyclic 

Left  rotor  lateral 
cyclic 

Elevator  incidence 

Lateral  cyclic  stick 

Rt  rotor  lateral 
cyclic 

Left  rotor  lateral 
cyclic 

Right  rotor  collec¬ 
tive 

Left  rotor  collec¬ 
tive 

Rt  rotor  F/A  cyclic 
Left  rotor  F/A 
cyclic 

Aileron 

Pedal 

Rudder  incidence 

Rt  rotor  collective 
Left  rotor  collec¬ 
tive 

Rt  rotor  F/A  cyclic 
Left  rotor  F/A 
cyclic 

Mast  tilt 

Collective  lower 
limit 

Collective-lateral 
cyclic  ratio 
Collective- rudder 
ratio 

TABLE  III.  SURFACE  TO  CONTROLS  ] 

Configuration 

Surface 

Control 

Single  Rotor 

Main  rotor  collective 

Collective  stick 

F/A  cyclic 

F/A  stick 

Lateral  cyclic 

Lateral  stick 

Tail  rotor  collective 

Pedal 

Rudder  incidence 

Pedal 

Aileron 

Lateral  stick 

Elevator  incidence 

F/A  stick 

Collective  stick 

Wing  incidence 

Collective  stick 

Jet  thrust 

Collective  stick 

Tandem 

Fwd  rotor  collective 

Collective  stick 

F/A  cyclic  stick 

Aft  rotor  collective 

Collective  stick 

F/A  cyclic  stick 

Fwd  rotor  F/A  cyclic 

F/A  cyclic  stick 

Aft  rotor  F/A  cyclic 

F/A  cyclic  stick 

Fwd  rotor  lateral  cyclic 

Lateral  cyclic  stick 
pedal 

Aft  rotor  lateral  cyclic 

Lateral  cyclic  stick 
pedal 

Rudder  incidence 

Pedal 

Wing  incidence 

Collective  stick 

Elevator  incidence 

Collective  stick 

F/A  cyclic  stick 

Jet  thrust 

Collective  stick 

Aileron 

Lateral  cyclic  stick 
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TABLE  III.  Continued 


Configuration 


Control 


Composite 


Rt  rotor  collective 


Collective  stick 
Lateral  cyclic  stick 
Pedal 
Mast  tilt 


Left  rotor  collective 


Collective  stick 
Lateral  cyclic  stick 
Pedal 
Mast  tilt 


Rt  rotor  F/A  cyclic 


F/A  cyclic  stick 
Lateral  cyclic  stick 
Pedal 


Left  rotor  F/A  cyclic 


F/A  cyclic  stick 
Lateral  cyclic  stick 
Pedal 


Rt  rotor  lateral  cyclic 


F/A  cyclic  stick 
Lateral  cyclic  stick 

Left  rotor  lateral  cyclic  F/A  cyclic  stick 

Lateral  cyclic  stick 

Rudder  incidence 


Pedal 


Elevator  incidence 


Collective  stick 
F/A  cyclic  stick 


Wing  incidence 
Aileron 


Collective  stick 


Lateral  cyclic  Stic! 
Collective  stick 


Jet  thrust 
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ERRATA 

The  following  corrections  are  applicable  to  AFFDL  TR  69-123,  A  Stability  and 
Control  Prediction  Method  for  Helicopters  and  Stoppable  Rotor  Aircraft,  Volume 
II  (User's  Manual),  February  1970. 

The  numbers  in  Figures  11  through  lU  which  are  associated  with  the  main  and 
tail  rotor  blade  twist  and  flapping  stop  should  be  negative.  The  correct  sign 
convention  is  illustrated  in  Figure  l8,  page  123.  To  correct  the  figures,  place 
a  negative  sign  in  front  of  the  "words"  as  indicated  below: 


Figure 

Page 

Group 

Line 

Word 

11 

112 

Main  Rotor 

1 

7 

11 

112 

Main  Rotor 

5,  6,  7 

1  through  7 

11 

112 

Tail  Rotor 

’  1 

7 

12 

lilt 

Main  Rotor 

1 

7 

12 

lilt 

Main  Rotor 

5.  6,  7 

1  through  7 

12 

lilt 

Tail  Rotor 

1 

7 

13 

116 

Main  Rotor 

1 

7 

13 

116 

Main  Rotor 

5,  6,  7 

1  through  7 

13 

116 

Tail  Rotor 

1 

7 

13 

116 

Tail  Rotor 

5.  6,  7 

1  through  7 

ll» 

118 

Main  Rotor 

1 

7 

lit 

118 

Main  Rotor 

5,  6,  7 

1  through  7 

lit 

118 

Tail  Rotor 

1 

7 

lU 

118 

Tail  Rotor 

5.  6,  7 

1  through  7 

The  computer  program  in  Volume  IV  of  AFFDL  TR  69-123  is  correct  as  printed. 
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